TOPICAL REVIEW 



The Gravitational Wave Signature 
of Core-Collapse Supernovae 

Christian D Ott 

E-mail: cottatapir.caltech.edu 

Theoretical Astrophysics, Mailcode 130-33, 

California Institute of Technology, Pasadena, California 91125, USA 
and 

Niels Bohr International Academy, Niels Bohr Institute, 
Blegdamsvej 17, 2100 Copenhagen, Denmark 

Abstract. We review the ensemble of anticipated gravitational-wave (GW) 
emission processes in stellar core collapse and postbounce core-collapse supernova 
evolution. We discuss recent progress in the modeling of these processes and 
summarize most recent GW signal estimates. In addition, we present new results 
on the GW emission from postbounce convective overturn and protoneutron star 
g-mode pulsations based on axisymmetric radiation-hydrodynamic calculations. 
Galactic core-collapse supernovae are very rare events, but within 3 — 5 Mpc from 
Earth, the rate jumps to 1 in ~ 2 years. Using the set of currently available 
theoretical gravitational waveforms, we compute upper-limit optimal signal-to- 
noise ratios based on current and advanced LIGO/GEO600/VIRGO noise curves 
for the recent SN 2008bk which exploded at ~3.9 Mpc. While initial LIGOs 
cannot detect GWs emitted by core-collapse events at such a distance, we find 
that advanced LIGO-class detectors could put significant upper limits on the 
GW emission strength for such events. We study the potential occurrence of 
the various GW emission processes in particular supernova explosion scenarios 
and argue that the GW signatures of neutrino-driven, magneto-rotational, and 
acoustically-driven core-collapse SNe may be mutually exclusive. We suggest 
that even initial LIGOs could distinguish these explosion mechanisms based on 
the detection (or non-detection) of GWs from a galactic core-collapse supernova. 



PACS numbers: 97.60.Bw, 97.60.Jd, 97.60.-s, 97.10. Kc, 04.30.Db, 04.40.Dg 

1. Introduction 

Ever since the very first experimental efforts to detect gravitational waves (GWs), core- 
collapse supernovae (SNe) have been considered as potential astrophysical emission 
sites. There are very strong indications from theory and observation that multi-D 
dynamics play a prominent and probably decisive role in core-collapse SNe (see, e.g., 
[1, 2]). GWs are emitted at lowest order by an accelerated mass-energy quadrupole 
moment. Hence, by their intrinsic multi-D nature, GWs, if detected from a core- 
collapse event, will very likely prove powerful messengers that can provide detailed 
and live dynamical information on the intricate multi-D dynamics occurring deep 
inside collapsing massive stars. 

Massive stars (8-10 Mq < M < 100 Mq at zero-age main sequence [ZAMS]) 
form electron-degenerate cores composed primarily of iron-group nuclei in the final 
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stages of their exoergic nuclear burning. Once such an iron core exceeds its effective 
Chandrasekhar mass (see, e.g., [3, 4]) it grows gravitationally unstable. Collapse 
ensues, leading to dynamical compression of the inner core material to nuclear 
densities. There, the nuclear equation of state (EOS) stiffens, resulting in the rebound 
of the inner core ( "core bounce" ) . A hydrodynamic shock wave is launched at the outer 
edge of the inner core and propagates outward in mass and radius, slamming into the 
still infalling outer core. Owing to the dissociation of heavy nuclei and to energy 
losses to neutrinos that stream away from the postshock region, the shock quickly 
loses energy, stalls and must be revived to plow through the stellar envelope, blow up 
the star, and produce a SN explosion, leaving behind a neutron star. Without shock 
revival, black-hole (BH) formation is inevitable and even with a successful explosion, 
a BH may still form via fall-back accretion. 

Iron core collapse is the most energetic process in the modern universe, liberating 
some 10^'^ erg ~ 100 B (Bethe) of gravitational energy. Most of this energy, ~ 99%, 
is emitted in neutrinos as the protoncutron star (PNS) contracts and cools over 
a timescale of 100 s. Only ^ 1% goes into the asymptotic energy of the SN 
explosion and becomes visible in the electromagnetic spectrum. The fundamental 
question that core-collapse SN theory has been facing for the past 45 years is 
how exactly the necessary fraction of gravitational energy is transferred to revive the 
shock and ultimately unbind the stellar envelope. Shock revival must occur sooner 
than 1-1.5 s after bounce (depending on progenitor star structure setting the rate of 
mass accretion) in order to produce a compact remnant that obeys observational and 
theoretical neutron star upper baryonic mass limits aroimd ~ 1.5-2.5 Mq (see [5] and 
references therein). 

The SN explosion mechanism may involve (a combination of) heating of the 
postshock region by neutrinos, multi-dimensional hydrodynamic instabilities of the 
accretion shock, in the postshock region, and in the PNS, rotation, PNS pulsations, 
magnetic fields, and nuclear burning (for a recent review, see [2], but also [6]). Three 
SN mechanisms are presently discussed in the literature. The neutrino mechanism has 
the longest pedigree [2, 4, 7], is based on postbounce neutrino energy deposition behind 
the stalled shock, and appears to require [8-10] convection and the standing-accretion- 
shock instability (SASI, see, e.g., [11] and references therein) to fimction in all but the 
very lowest-mass massive stars which may explode even in spherical symmetry [12, 13]. 
Recent detailed 2D neutrino-radiation hydrodynamics simulations by the Garching 
group produced explosions in particular 11.2-A/q and 15-A/0 progenitor models 
[8, 9]. However, it is not yet clear how the neutrino mechanism's efficacy varies 
with progenitor ZAMS mass and structure and what its detailed dependence on the 
high-density nuclear EOS may be. 

The magneto-rotational (or MHD) mechanism, probably operating only in the 
context of rapid progenitor rotation, depends on magnetic-field amplification during 
collapse and at postbounce times. It leads to explosions that develop in jet-like fashion 
along the axis of rotation [14-18] and may reach hypernova energies of ~ 10 B [17]. 
The MHD mechanism may also be relevant in the context of long-soft gamma-ray 
bursts (GRBs, see, e.g., [19]) and could be a precursor, setting the stage for a later 
GRB (e.g., [17], but see [20]). 

The acoustic mechanism for core-collapse SNe, as recently proposed by 
Burrows et al. [6, 21-23], requires the excitation of large-amplitude PNS pulsations 
(primarily g-modes) by turbulence and SASI-modulatcd accretion downstreams. 
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These pulsations damp by the emission of strong sound waves that steepen to shocks 
and deposit energy in the postshock region, eventually leading to late explosions at 
> 1 s after bounce. This mechanism appears to be sufficiently robust to blow up even 
the most massive and extended progenitors, but has so far not been confirmed by 
other groups (see also [24, 25]). 

Constraining the core-collapse SN mechanism via astronomical observations is 
difficult. The intricate pre-explosion dynamics of the SN core deep inside the 
supergiant presupernova star are inaccessible by the traditional means of astronomy. 
Theoretical models of the SN mechanism can currently be tested via secondary 
observables only, including the asymptotic explosion energy, ejecta morphology, 
nucleosynthesis products, compact remnant mass and proper motion, and pulsar 
spin/magnetic fields. 

GWs and neutrinos are the only messengers with the potential of delivering first- 
hand information on the physical processes leading to explosion: Both are emitted deep 
inside the SN core and travel to observers on Earth practically without interaction 
with intervening material. 

A small number of neutrinos were detected from SN 1987A in the Large 
Magellanic (distance D « 50 kpc; see, e.g., [4] and references therein). GWs have 
not yet been observed directly, but the advent of GW astronomy has begun. An 
international network of broad-band light-interferometric GW observatories is active, 
encompassing the US LIGOs [26], the British-German GEO600 [27], the French- Italian 
VIRGO [28] and the Japanese TAMA 300 [29]. The three LIGO interferometers have 
recently reached their design sensitivities, and, in their S5 science run, have taken a 
year worth of data, partly in coincidence with GEO600 and VIRGO. In addition, there 
are a number of active resonant bar/sphere GW detectors in operation, including the 
four bar detectors of the International Gravitational Event Collaboration (IGEC-2), 
ALLEGRO, AURIGA, EXPLORER, and NAUTILUS [30], and the resonant spheres 
MiniGrail [31] and Schenberg [32]. The current status of ground-based GW detection 
was recently summarized by Whitcomb [33] . 

GWs from astrophysical sources are weak and notoriously difficult to detect (e.g., 
[34]). Hence, in order to disentangle an astrophysical GW signal from the mostly 
overwhelming detector noise, GW astronomy does not only require sensitive detectors, 
but also extensive processing and analysis of the detector output on the basis of reliable 
theoretical estimates for the GW signals presently expected from astrophysical sources. 
The latter must, in most cases, be obtained via detailed numerical modeling of the 
dynamics responsible for the GW emission in a given source. 

In iron core collapse and postbounce SN evolution, the emission of GWs is 
expected primarily from rotating collapse and bounce, nonaxisymmetric rotational 
instabilities, postbounce convective overturn/SASI, and PNS pulsations. In addition, 
anisotropic neutrino emission, global precollapse asymmetries in the iron core and 
surrounding burning shells, aspherical mass ejection, magnetic stresses, and the late- 
time formation of a black hole may contribute to the overall GW signature. 

The aim of this topical review is to summarize the recent significant progress 
in the modeling of the various GW emission processes in core-collapse SNe with a 
focus on the early postbounce, pre-explosion SN evolution up to ~ 1 — 2 seconds 
after bounce. We do not cover the GW emission from the collapse of supermassive 
primordial or very massive population III stars, dynamical fission processes, late-time 
fall-back accretion on black holes, nuclear phase-transitions in PNSs, or from late-time 
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postbounce secular instabilities such as y'-modcs. Other reviews of the GW signature 
of core-collapse SNe are those of Kotake et al. [35] and Fryer and New [36] . 

In section 2, we provide a concise historical overview on early work and go on to 
discuss computational core-collapse SN modeling and the various ways in which GR 
and GW extraction arc treated in section 3. Section 4 covers the most extensively 
modeled GW emission mechanism, rotating core collapse and bounce. The potential 
for and the GW emission from nonaxisymmetric rotational instabilities is the topic 
of section 5. Section 6 is devoted to the emission of GWs from convective overturn 
and SASI, while we discuss the GW signal from PNS core pulsations in section 7. 
To both of these sections we add new, previously unpublished results obtained via 2D 
Newtonian radiation- hydrodynamics calculations. In section 8, we discuss the emission 
of GWs from anisotropic neutrino radiation fields and in section 9 we summarize the 
GW signals associated with rapid aspherical outflows, precollapse global asymmetries, 
strong magnetic stresses, and PNS collapse to a black hole. 

While the SN rate in the Milky Way and the local group of galaxies is rather 
low and probably less than 1 SN per two decades (e.g.. [37]), there may be 1 SN 
occurring about every other year between 3 — 5 Mpc from Earth [38] . The recent SN 
2008bk, which exploded roughly 3.9 Mpc away, is an example SN from this region of 
space. Thus, in section 10, we present optimal single-detector matched-filtering signal- 
to-noise ratios for LIGO/GEO600/ VIRGO and advanced LIGOs for a subset of the 
gravitational waveforms reviewed in this article and with an assumed source distance 
corresponding to that of SN 2008bk. We find that initial LIGO-class detectors had no 
chance of detecting GWs from SN 2008bk. Advanced LIGOs, however, could put some 
constraints on the GW emission strength, but still would probably not allow detailed 
GW observations. We wrap up our review in section 11 with a critical summary of 
the subject matter and discuss in which way the various GW emission processes can 
be linked to particular SN explosion mechanisms. We argue that the GW signatures 
of the neutrino. MHD. and acoustic SN mechanisms may be mutually exclusive and 
that the mere detection, or, in fact non-detection, of GWs from a nearby core-collapse 
SN can constrain significantly the core-collapse SN explosion mechanism. 

In this review article, all values of the dimensionless GW signal amplitudes h+ 
and hx are given for optimal source-observer orientation and a source distance of 
10 kpc is typically assumed. In most figures showing GW signals, the waveforms are 
plotted as /i+,x D, rescaled by distance D and in units of centimeters. Summaries of 
GW extraction methods can be found in [34, 39-42] and a comparison of Newtonian 
and general relativistic methods were presented in [43, 44]. 

2. A Short Overview on Early Work 

Although gravitational waves were proposed by Einstein himself in 1918 [45], there 
was general disagreement on the waves' reality and their ability to carry off energy 
from their source until a 1957 gedankenexperiment by Bondi [46] restored trust in the 
fundamental theory of GWs (see, e.g., [34] and references therein for more details). 
This lead to a renaissance of GW physics in the early to mid 1960s, resulting in much 
theoretical work on astrophysical sites of GW emission (see, e.g., [34. 47]) and to the 
first experimental searches for GWs led by Weber with resonant bar detectors [33, 48]. 

GW emission in the core-collapse SN context first appeared in a 1966 paper 
by Weber [49] in which he referred to unpublished work from 1962 by Dyson [50] 
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on astrophysical sources of GWs, including binary neutron star systems and non- 
spherical stellar collapse^. Also in 1966, Wheeler, based on unpublished work by 
Zee and himself, published first quantitative estimates of the emission of a GW 
burst by quadrupole oscillations of a PNS associated with core bounce [52]. In the 
following years, the theory of NS pulsations and their GW emission characteristics 
was worked out by Thorne, Bardeen, Meltzer, and collaborators [53-59], leading to 
many subsequent studies and opening up a rich independent branch of research (see, 
e.g., [60, 61] for reviews). 

In their 1971 review [51], Ruffini and Wheeler provided a first comprehensive 
summary of scenarios for GW emission in the stellar collapse context and provided 
quantitative estimates for the GW signal from rotating core collapse and bounce, (early 
postbouncc) neutron star pulsations, and nonaxisymmctric neutron star deformations. 

A number of subsequent studies in the 1970s and early 1980s focussed on 
aspherical collapse and relied upon semi-analytical descriptions of collapsing and/or 
oscillating homogeneous spheroids and ellipsoids. Key studies were those carried out 
by Thuan and Ostriker [62], Novikov [63], Shapiro [64], Saenz and Shapiro [65- 
67], Epstein [68], Moncrief [69], and Detwciler and Lindblom [70]. Most of 
these early and ground-breaking investigations severely overestimated the strength 
(amplitude/energy) of the GW burst from core collapse and bounce and failed to 
capture the qualitative features of the core bounce GW burst observed later in 
numerical simulations. Some of the former and other authors using similar techniques 
also studied the GW emission from nonaxisymmctric rotational deformations of 
rapidly rotating postbouncc neutron stars [67, 71-75]. 

In 1979, Turner and Wagoner [76] followed a different approach and used 
perturbation theory to calculate the GW emission produced by rotationally-induced 
perturbations of spherically symmetric numerical iron core collapse models of Van 
Riper [77] and Wilson [78]. Seidcl and collaborators followed along those lines and 
perturbatively analyzed detailed spherically-symmetric general-relativistic (GR) core 
collapse simulations in the 1980s [79, 80]. Although the perturbative approach is 
in principle appropriate only for small deviations from sphericity and, hence, for very 
slowly rotating stellar cores, the GW signals from core bounce predicted by Seidel et al. 
is in rough qualitative agreement with the signals obtained from more recent multi-D 
simulations of even rapidly rotating stellar cores. 

Epstein [81] and independently Turner [82] considered in their 1978 studies the 
GW emission from anisotropic neutrino emission in stellar core collapse and developed 
the linear-theory formalism that has since been applied to extract GWs from aspherical 
neutrino radiation fields in multi-D radiation-hydrodynamic simulations. 

The early epoch of work on the GW emission in stellar core collapse was 
dominated by analytic, semi-analytic, and perturbative studies, largely because of 
a lack of computers sufficiently powerful to carry out multi-D numerical non-linear 
hydrodynamics simulations. The end of this epoch was marked by the emergence of 
the first supercomputers and their general availability to the astrophysics community 
beginning in the late 1970s and early 1980s. 

X Upon being asked by the author to comment on his early work on GW emission from stellar 
collapse, Dyson wrote: / have no idea who first calculated the emission of gravitational radiation 
from a collapsing star with rotation. Ruffini and Wheeler [51] may have been the first. It was 
certainly obvious to everyone who thought about it that a collapsing star with rotation would give 
rise to a strong pulse of gravitational waves. I make no claim to have thought of this first. F. Dyson, 
priv. comm., Oct. 2006. 
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3. Modeling Core- Collapse Supernovas 

The Einstein equations as well as the equations of radiation-magnctohydrodynamics 
are families of coupled non-linear partial differential equations for which solutions can 
be found analytically only in very few limiting cases, leaving numerical modeling as 
the only viable route to their solution in complex scenarios such as stellar collapse 
and core-collapse SNe. Hence, for an accurated understanding of the core-collapse SN 
phenomenon and the GW signals emitted by the involved multi-D dynamics, detailed 
numerical models are required. 

Modeling stellar core collapse and the postbounce evolution of the SN core is a 
multi-scale multi-physics problem that involves lengthscales from the extended pre-SN 
stellar core (thousands of kilometers) down to small-scale turbulence in the postbounce 
flow (on the order of meters) and timescales from < 10~^ s (the typical computational 
timestep) to up to 1 — 2 s (the time for the development of a full explosion or for BH 
formation to occur). 

An ideal and complete ("realistic") model should resolve all lengthscales, should 
be free of artificially imposed symmetries (i.e., be 3D) and should fully include general 
relativity, GR MHD, multi-D GR Boltzmann neutrino transport, and a microphysical 
treatment of EOS and nuclear and neutrino interactions. In addition, since the high- 
density EOS, the pre-SN stellar structure and angular momentum configuration are 
not well constrained; extensive parameter studies are necessary to plow the parameter 
space of possible EOS and pre-SN configurations. Present-day 2D core-collapse SN 
simulations, e.g., those carried out by [9, 17, 22, 83], scale to ^ 48-128 compute 
cores and typically require multiple months to complete on current supercomputers. 
A single 3D full GR radiation- (magneto)hydrodynamics simulation will have an at 
least a 100 times (and more realistically, 1000 to 10000 times) greater computational 
complexity and will need to scale to thousands of compute cores on next-generation 
supercomputers to complete within reasonable time§. 

None of the currently published numerical models of core collapse and postbounce 
SN evolution live up to the above standards and all studies employ approximations 
for multiple of the physical ingredients listed. A detailed discussion of the various 
approximations cannot be provided here, but since the focus of this review is on GWs, 
a consequence of GR and. in principle, requiring a full GR treatment, we present in 
the following the different ways in which GR gravity, relativistic dynamics, and GW 
extraction from matter and spacetime dynamics are presently handled in multi-D 
core-collapse simulations. 

(i) Newtonian gravity and dynamics. Many recent multi-D simulations 
(e.g. ,[17, 21-23, 41, 83, 85-90]) employ a solution of the Newtonian Poisson 
equation for computing the gravitational acceleration terms and treat (magneto)- 
hydrodynamics and radiation transport (if included) in Newtonian fashion. GWs 
are extracted from the fluid motions via the slow-motion, weak-field quadrupole 
formalism (e.g., [34, 91, 92]). 

(ii) Approximate GR gravity. Motivated by the work of the Garching SN 
group [93-95], a set of recent core-collapse SN simulations [8, 9, 96-98] 
take GR into account via replacing the spherical component of the multipole 
decomposition of the Newtonian gravitational potential with an effective GR 

§ See [84] for a more quantitative discussion of the computational complexity of the core-collapse SN 
problem. 
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potential modeled after the Tolman-Oppenhcimer-Volkoff (TOV) potential (e.g., 
[99]). The hydrodynamics, however, are treated in standard Newtonian fashion. 
Redshift and time-dilation effects are taken into account in the radiation transport 
sector in the simulations of the Garching SN group [8, 95]. GWs are extracted 
via the slow-motion, weak-field quadrupole formalism [34, 91, 92]. 

(iii) Conformally-fiat (CFC) GR. In the conformally-flat (or conformal-flatness 
condition [CFC]) approximation to GR introduced by Isenbcrg [100] and first 
used by Wilson ct al. [101], the general 3-mctric of the 3 + 1 decomposition (e.g., 
[102]) is replaced by the flat-space Minkowski metric scaled with a conformal- 
factor 0^. The CFC approximation is exact in spherical symmetry. In multi-D, 
a CFC spacetime behaves as an approximation of GR at first post-Newtonian 
order [103] and may be regarded as GR minus the dynamical degrees of freedom 
of the gravitational field that correspond to the GW content at infinity. However, 
even stationary spacetimes without GWs can be non-CFC, e.g., rotating NSs 
or Kerr BHs (e.g., [104, 105]). An extension of the CFC approach with terms 
of second post-Newtonian order, CFC-I-, was introduced by Ccrda-Duran and 
collaborators [106]. The CFC approximation permits the use of the full GR 
(M)IID equations (e.g., [107]) and is employed in simulations of axisymmctric 
rotating core collapse by Dimmelmcicr and collaborators [108--114]. At least 
for this scenario it has been shown to be an excellent approximation to full GR 
[106, 115]. Since a CFC spacetime does not contain GWs, they must be extracted 
from the fluid dynamics via variants of the slow-motion, weak-field quadrupole 
formalism, making use of the conserved variables of GR (M)IID (e.g., [108, 112]). 

(iv) Full GR. Multi-D simulations in full 2-1-1 and 3 + 1 GR evolve the Einstein 
equations without approximations and with appropriate choices of temporal and 
spatial gauge (e.g., [102]). Full GR simulations of rotating core coUapse and 
early postbounce evolution were carried out by Shibata et al. [116-118] and by 
Ott et al. [42, 115, 119] who employed the BSSN spacetime evolution system 
[102, 120, 121]. In principle, GWs can be extracted directly from the spacetime 
in regions sufficiently far away from the emission region (e.g., [122] and references 
therein). Due to the relative weakness of the GWs emitted in core collapse, the 
direct extraction has proven to be numerically difficult and most studies resort to 
the same slow- motion, weak- field quadrupole formalism used in the CFC context ||. 

Studies addressing the GW signature of stellar core collapse and core-collapse SNe 
historically followed two general approaches. 

One approach focussed on the detailed treatment of EOS, microphysics and 
radiation transport as necessary for physically accurate models of the collapse and 
postbounce phase. Studies following this approach were either Newtonian or employed 
approximate GR, neglecting the additional complication of CFC or full GR (e.g., 
[23, 39, 86, 87, 97, 123, 124]). Due to the large computational complexity (dominated 
by neutrino transport) of such detailed simulations, these studies were limited to small 
model sets and, hence, were unable to cover the parameter space of possible initial 
configurations and EOS. 

Studies following the second approach were primarily concerned with the GW 
signal of rotating core collapse and bounce. Their emphasis was on performing large 

II See [43, 44] for comparisons of gauge-invariant extraction techniques with variants of the quadrupole 
formalism. 
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Figure 1. Axisymmetric GW burst signal {h-^ D in units of cm, where D is the 
distance of the source) as a function of time after core bounce for models A1B3G3 
(type I), A2B4G1 (type II), and A1B3G5 (type III) of the 2002 GR study using 
polytropic initial models and a simple analytic hybrid polytropic/ideal-fluid EOS 
carried out by Dimmelmeier et al. [112]. The waveforms were obtained from [126]. 



parameter studies of precollapse rotational configurations in Newtonian gravity, CFC, 
or full GR and, in some cases, included magnetic fields. Their schemes typically 
did not take into account a microphysical EOS, detailed neutrino transport and 
microphysics, and presupernova models from stellar evolutionary calculations (see, 
e.g., [40, 41, 90, 112, 116, 125], but note that [41] employed a microphysical EOS, but 
neglected deleptonization) . 

Only recently, these two approaches have begun to come together in the 
first extensive CFC and full GR parameter studies of rotating core collapse with 
a microphysicsl EOS and an approximate deleptonization scheme carried out by 
Dimmelmeier et al. [108, 109] and Ott et al. [42, 115, 119]. The results of these 
studies will be discussed extensively in the following section. 

4. Rotating Core Collapse and Bounce 

Rapid precollapse rotation, in combination with angular momentum conservation 
during collapse, leads to significant asphericity in the form of a GW emitting, rapidly 
time-varying £ = 2 oblate (quadrupole) deformation of the collapsing and bouncing 
core. 

Rotating collapse and core bounce is the most extensively studied GW emission 
process in the massive star collapse context. In 1982, Miiller published in [127] 
the first GW signals from rotating core collapse and bounce that were based on 
the axisymmetric (2D) Newtonian simulations of Miiller and Hillebrandt [128]. A 
large number of 2D studies followed with varying degrees of microphysical detail, 
inclusion of GR, and precollapse model sets (see, e.g., [39-41, 86, 112, 116, 129] and 
references therein). These computational parameter studies demonstrated that the 
general analytic picture of stellar core collapse derived by Goldreich & Weber [130] 
and Yahil [131] for spherically-symmetric collapse also holds to good approximation 
for rotating cores: From the beginning of collapse, the collapsing core separates into a 
subsonically homologously (f cx r) contracting inner core and a supersonically infalling 
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outer core. The mass of the inner core at core bounce sets the mass cut for the matter 
that is dynamically relevant in bounce (see, e.g., [4]) and responsible for the GW 
burst. It also determines the initial size of the PNS. 

Furthermore, these studies identified at least three GW signal "types" that can 
be associated with distinctly different types of collapse and bounce dynamics (figure 1 
displays representative examples). Type I models undergo core bounce governed 
by the stiffening of the nuclear EOS at nuclear density and "ring down" quickly 
into postbounce equilibrium. Their waveforms exhibit one pronounced large spike 
at bounce and then show a gradually damped ring down. Type II models, on the 
other hand, are affected significantly by rotation and undergo core bounce dominated 
by centrifugal forces at densities below nuclear. Their dynamics exhibits multiple 
slow harmonic-oscillator- like damped bounce-re-expansion-collapse cycles ( "multiple 
bounces" ) , which is reflected in the waveform by distinct signal peaks associated with 
every bounce. It is interesting to note that type-II models are related to fizzlers, 
proposed (e.g., [132-134]) collapse events of self-gravitating fluid bodies that become 
temporally or permanently stabilized by centrifugal forces before core bounce or black- 
hole formation. 

Type III models are characterized by fast collapse (owing to a very soft subnuclear 
EOS or very efficient electron capture), extremely small masses of the homologously 
collapsing inner core, and low-amplitude GW emission and a subdominant negative 
spike in the waveform associated with bounce. 

GWs from collapse with magnetic fields were studied by [87, 90, 96, 117, 125] who 
found an additional dynamics/signal type IV which occurs only in the case of very 
strong precollapse core magnetization {B > 10^^ G). Such strong precollapse magnetic 
fields are unlikely to occur in nature and the aforementioned studies showed that 
weaker magnetic fields have little dynamical consequence during collapse and bounce 
(see also the discussion in [17]). However, MHD-drivcn jet-like bipolar outflows and 
time-changing magnetic fields themselves can lead to the emission of GWs (see the 
discussion in section 9) 

Recently, Ott et al. [115, 119] and Dimmelmeier et al. [108, 109] presented the first 
3D full GR and 2D CFG GR simulations of rotating stellar collapse that incorporate 
GR as well as a microphysical finite-temperature nuclear EOS and an approximate 
scheme by Liebendorfer (first presented in [139]) to account for electron capture and 
neutrino losses during collapse. These calculations for the first time included all the 
known most relevant physics in the collapse and bounce phase and, in the case of 
[108], also considered two different nuclear EOS, the Lattimcr-Swesty EOS [140] and 
the Shen EOS [135, 136]. 

The results of Dimmelmeier et al. [108] show that the combined effects of 
deleptonization and GR decrease the mass and extent of the homologous inner core 
and reduce significantly the relevance of centrifugal support for a large set of progenitor 
models and a wide range of precollapse rotational configurations. This set encompasses 
and goes beyond what is deemed realistic in the context of collapsing massive stars 
(see, e.g., [141, 142]). In particular, they find that the GW signal of rotating collapse 
and core bounce is of generic Type-I shape. Type-II dynamics with the associated 
"multiple bounce" GW signals do not obtain when GR and deleptonization are 
included. Type-III dynamics are also absent, since they require very efficient electron 
capture and resulting very small masses of the inner core < 0.3 Mq that do not occur 
in the iron core collapse context where the smallest inner core mass is > 0.45 Mq [108]. 
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Figure 2. GW signals D in units of cm, where D is the distance of the source) 
for a few examples from the 2D GR model set of Dimmelmeier et al. [108]. The 
models shown here were computed with the Shen EOS [135, 136] and employ ID 
presupernova models of [137], spanning the progenitor mass range from 11.2 Mq 
(sll) to 40 Mq (s40). The models were set up with precollapse central angular 
velocities ; from ~ 1.5 rads~^ to ~ 11 rads~^. For details of the rotational 
setup, see [108] . Model E20A uses a 20-A'/q presupernova model that was evolved 
by [138] with a ID prescription for rotation. Note the generic shape of the 
waveforms, exhibiting one pronounced spike at core bounce and a subsequent ring 
down. Very rapid precollapse rotation {Qc.i ^ 6 rads~^; models s20A3O12 and 
s40A2O13 in this plot) results in a significant slow-down of core bounce, leading 
to a lower-amplitude and lower-frequency GW burst. The GW signal data are 
available for download from [126]. 



Table 1 . Summary of the G W signal characteristics of rotating iron core collapse 
and core bounce based on the waveforms of Dimmelmeier et al. [108]. All models 
exhibit type-I dynamics and waveform morphology and can be organized into three 
distinct groups based primarily on their precollapse central angular velocity 
l^max] is the maximum gravitational wave strain amplitude (scaled to 10 kpc) 
at bounce, -Egw is the energy radiated away in gravitational waves, /peak is the 
frequency at which the GW energy spectrum dEQ^fj/df peaks, and A/50 is the 
frequency interval centered about /peak that contains 50% of Eq^/. Note that 
/max used by Dimmelmeier et al. is the peak of the GW signal spectrum and 
not the peak of dEQ^fj / df . Also note that for the slowly rotating group prompt 
postbounce convective overturn contributes significantly to the overall GW signal. 
The convective contribution was removed from the waveforms before analysis, 
since the delcptonization scheme employed by Dimmelmeier et al. and Ott et al. 
is ineffective at postbounce times and is likely to overestimate the strength and 
duration of prompt convection after bounce [108]. 



Group 
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(rad s"-') 


(10^21 at 10 kpc) 
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(Hz) 


(Hz) 


1 


< 1-1.5 


£0.5 


<o.i 


~ 700-800 


~ 400 


2 


1-2 to 6-13 


0.5 to 10 


0.1 to 5 


~ 400-800 
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3 > 6-13 3.5 to 7.5 0.07 to 0.5 70 to 200 80 to 250 
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Furthermore. Dimmelmeicr et al. [108] demonstrate that the rotating collapse 
and bounce dynamics and the resulting GW burst signal depend primarily on the 
precollapse central angular velocity f2c,i £^nd secondarily on the progenitor mass, which 
influences precollapse iron core entropy and mass, which, in turn affects the mass of 
the inner core at bounce (see, e.g., [3, 143]). The degree of differential rotation and 
the stiffness of the particular choice of nuclear EOS appear to have little influence on 
the general characteristics of the GW signal. 

Figure 2 shows example waveforms of a subset of the models considered by 
Dimmelmeier et al. [108] whose waveform data are available from [126]. Although 
all models exhibit generic type-I dynamics and signal morphology, the dynamics and 
GW signals can be ordered into three groups (for more details, see [108] and the 
summary in table 1): 

(1) Slowly rotating iron cores with ^lc,i ^ 1-1.5 rads~^ undergo core bounce 
dominated by nuclear pressure effects, develop only a small quadrupole 
deformation and yield small peak GW amplitudes [/imaxl at core bounce that are 
generally below ~ 5 x 10~^^ for a galactic core-collapse event at 10 kpc distance. 
All slowly rotating models exhibit strong prompt postbounce convective overturn 
in the postshock region owing to the negative entropy gradient left behind by 
the stalling the shock. Because of the lack of postbounce neutrino transport, 
the strength and duration of the prompt convection is probably overestimated in 
these models [108] (see also section 6.1 in this article). 

(2) Moderately rapidly to rapidly rotating iron cores with 1-2 rads~^ < ric,i ^ 
6-13 rads~^ still experience pressure-dominated bounce, but develop larger 
quadrupole deformations and have a rotationally increased mass of the inner core 
at bounce. This results in sizeable peak GW amplitudes and energy emissions 
(see table 1). Prompt postbounce convection is suppressed by positive specific 
angular momentum gradients in the postshock region and does not contribute 
to the GW signal (see, e.g., [141]). The upper end of the range in Oc,i for this 
group is dependent on progenitor characteristics. Progenitors with massive (and 
higher-entropy) iron cores tend to transition to the next group at lower ilc^i- 
Dimmelmeier et al. find approximate transition fic,i of > 13, ~ 9, ^ 11, and 
^ 7 rad s^^ for their 11.2, 15, 20, and 40 M© models, respectively. 

(3) Very rapid precollapse rotation (ilc.i ^ 6-13 rads~^; the actual value depending 
on the progenitor model) results in slow core bounce governed exclusively by 
centrifugal forces at significantly subnuclear densities. This results in a decrease 
of |/imax| and a shift of the peak of the GW energy spectrum to frequencies below 
~ 400 Hz. 

Figure 3 contrasts the strengths of the GW bursts of the models of 
Dimmelmeier et al. [108] with initial LIGO sensitivity [26]. Plotted is the initial 
LIGO dimensionless rms strain sensitivity ft-rms = a/ fS{f) as a function of frequency 
/, where S{f) is the noise spectral density in units of (Hz)^^/^. Each dot represents a 
particular model and is plotted at that model's detector-dependent characteristic GW 
frequency fc and characteristic strain he, computed by Dimmelmeier et al. according 
to the prescription given in [34]. The models cluster in the hc-fc plane according to 
their membership in one of the three groups discussed above. Slowly rotating models 
(group 1) have small he and moderately high fc. With increasing ili^c, the model 
move upward and slightly to the right (higher he and slightly higher fc). Moderately 
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Figure 3. Location of the GW signals of Dimmelmeier et al. [108] in the hc-fc 
plane (using the detector-dependent definitions of he and fc given in [34]) relative 
to the /irms = V fS{f} sensitivity of initial LIGO and for a source distance of 
10 kpc. The three groups of signals discussed in the text and summarized in 
table 1 are marked. This plot is a similar to figure 16 of [108] and based on freely 
available data [126]. 

rapid models (group 2) reach high he and cluster in fc between 600 and 800 Hz. The 
most rapidly rotating models of that group move somewhat to the left in f^, while 
very rapidly rotating and ccntrifugally bouncing models (group 3) cluster at low 
and lower he, but still exhibit a high signal to noise ratio (SNR) with initial LIGO 
and at a distance of 10 kpc. Hence, according to figure 3, GWs from moderately 
rapidly to rapidly rotating core collapse and bounce within the Milky Way should be 
detectable by current LIGO-class detectors. In addition, it may be possible, based on 
a determination of the dominant emission frequency (ideally in combination with a 
determination of the progenitor mass by other astronomical means), to constrain the 
precollapse iron core rotation. 

To conclude this section, we point out that according to the most recent 
stellar evolutionary models and pulsar birth spin estimates [141, 142], garden- variety 
massive stars at solar metallicity are probably spinning with central angular velocities 
significantly below 1 rads~^ (group 1). GW bursts from core bounce of such models 
are unlikely to be detectable by current detectors and may be marginally detectable 
by advanced LIGO. More rapid rotation may be relevant in the context of gamma-ray 
burst progenitors (see, e.g., [19]) which could make up perhaps ^ 1% of the single 
massive star population, but may occur primarily at low metallicitics [144]. 

5. Rotational Nonaxisymmetric Instabilities 

From the classical theory of Newtonian MacLaurin spheroids (see, e.g., [145, 146]) one 
finds that rotating axisymmetric fluid bodies become unstable to nonaxisymmetric 
deformations due to a dynamical instability at rotation rates P = T/\W\ > 0.27 = 
Pdyn, where T is the rotational kinetic and \ W\ the gravitational energy of the spheroid. 
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A secular (dissipation driven^) instability may set in at /3 = > 0.14 = /3scc, 

though, due to its secular nature, has longer growth times than the dynamical 
instability. 

Stars that go unstable to the classical MacLaurin-typc dynamical or secular 
instability develop global azimuthal (nonaxisymmetric) structure that, at least in 
the linear regime, can be characterized in terms of modes m with spatial structure 
proportional to exTp{imip), where ip is the azimuthal angle. In most cases the 
bar-like m ~ 2 mode is dominant and one frequently speaks of a "barmode 
instability". Direct numerical simulations have demonstrated (e.g., [147-149]) that 
Pdyn holds approximately even in the case of differentially rotating compressible fluid 
configurations in GR. For the secular instability, perturbative studies suggest an onset 
at significantly lower /3 in GR (see, e.g., [61, 146] for discussions). However, GR 
hydrodynamic studies of the secular instability remain yet to be carried out. 

A spinning bar is a simple system that emits GWs at twice its rotational frequency 
(due to its TT-symmetry) with amplitudes |/ibar| oc MR^il^ /D, where M is the 
bar's mass, 2R its length, and fl ~ 211 f its angular velocity. Using the Newtonian 
quadrupole approximation, one can derive an estimate for the GW amplitude [88] , 

|(.U»«xio-'(jij)(^)% 

C^^-V— If— 

VlOkpcy V0-7MQyV12km 

where e is the ellipticity of the bar. The scaling of (1) is set up to reflect a PNS core 
of 0.7 Mq with a 2 ms period and a radius of 12 km that is deformed only moderately 
(e = 0.1) and located at a distance of 10 kpc. 

5.1. Instability at High- (3 

The literature on nonaxisymmetric instabilities in rotating fluid bodies is extensive and 
cannot be reviewed here in the necessary detail. Relatively recent reviews and relevant 
references can be found in [61, 146, 148]. In the context of massive star collapse and 
PNSs, the first Newtonian 3D simulations were carried out by Rampp et al. [150] and 
Brown [151] who followed the postbounce development of nonaxisymmetric dynamics 
in extremely differentially and rapidly rotating simplified stellar models that reached 
up to /3 « 0.35 at bounce in 2D and were mapped to 3D shortly before or after bounce. 
With a similar approach, though in full GR, Shibata & Sekiguchi [118] studied the 
development of nonaxisymmetric dynamics in models that reached (3 near /3dyn and 
reported the development of m = 2 and to = 1 dynamics in their models. 

Ott et al. [115, 119, 141] and very recently Dimmelmeier et al. [108] studied 
via their 3D full GR and 2D CFC GR core-collapse simulations the prospects for 
nonaxisymmetric rotational instabilities in massive star collapse and postbounce 
evolution. Their models included microphysical details and they considered a large 
number of precollapse rotational configurations that covered a parameter space, 
encompassing and going beyond current predictions from stellar evolution theory and 
constraints from pulsar birth spin estimates [141, 142, 144]. 

Ott et al. and Dimmelmeier et al. found that even their most extreme models 
do not reach values of /3 close to /3dyn during collapse, bounce, and during early 

X Possible mechanisms are physical viscosity and gravitational-wave reaction. See, e.g., [61, 146]. 
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Figure 4. Left; Colormap depicting the specific entropy distribution in tiie 
equatorial plane of the 20-Mo model s20A2B4 of Ott et al. [42, 115, 119] at 
90 ms after core bounce. Red and yellow regions of the PNS core have low 
entropy (~ 1-3 kfl/baryon) while dark blue and black symbolize high entropy 
> 6 ke/baryon. The nonaxisymmetric structures are of primarily m = {1,2,3} 
nature with radial variations in dominance. Right: GW polarizations h+ (top 
panel) and hx (bottom panel) multiplied by distance D as seen by observers 
along the equator (black) and along the pole (red) in the same model. Significant 
nonaxisymmetric dynamics with m = 2 components develop around 30—40 ms 
after core bounce at a rotation rate /3 ~ 0.13. Note that the GW burst signal from 
core bounce is purely axisymmetric, since an axisymmetric system has vanishing 
hx and vanishing GW emission along the axis of symmetry. The waveforms are 
available from [155]. 



postbouncc times. Their work confirmed results obtained in studies considering 
evolutionary sequences of rotating PNSs [152] and 2D Newtonian simulations [41, 141]. 

On the basis of the current understanding, it appears rather unlikely that a PNS 
in nature develops a high-/? rotational instability before, at, or early after bounce. 
As the PNS contracts and cools on a timescale of seconds to minutes and if angular 
momentum is conserved and not redistributed by other means (see, e.g., [17, 117, 141]), 
it spins up and may reach /3dyn- It is, however, more likely that the secular instability 
driven by viscous dissipation or GW back reaction, which in PNSs has a growth 
timescale of the order of or larger than ~ 1 s [153, 154], may set in first. 

5.2. Instability at Low-(3 

Via 3D full GR simulations that were carried out self-consistcntly from the onset 
of collapse, Ott et al. [42, 115, 119] confirmed that stellar cores with precollapse 
parameters in the investigated parameter space stay axisymmetric throughout 
collapse, bounce and the very early postbounce phase. Furthermore, they found 
the development of nonaxisymmetric dynamics at postbounce times > 20 ms in a 
subset of their models. The 3D dynamics were of predominantly m = {1, 2, 3} spatial 
character (see the left panel of figure 4) and were followed to late times in only two 
models, one having a postbounce (3 of 0.09, the other one spinning with f3 « 0.13. Very 
recently, Scheidegger et al. [96] , who performed a 3D Newtonian study with a monopole 
gravitational potential with relativistic corrections, found a similar instability in a 
model with a smaller postbounce value of P around 0.04-0.05. 
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The instability at low-/? seen in simulations may be related to a class of co-rotation 
dynamical shear instabilities that operates on the shear-energy stored in differential 
rotation and that powers azimuthal modes via resonant coupling with the background 
fluid. This occurs at co-rotation points where the mode pattern speed coincides with 
the fluid angular velocity [156, 157]. Such a low-/3 instability was first observed in a 
compact-star context by Centrella et al. [158] in 2001 and multiple subsequent studies 
confirmed its presence in differentially rotating equilibrium stellar models below the 
classical high-/3 instability limits (e.g., [159-162]). In 2005, Ott et al. [163] discovered 
a low-/3 instability in a 3D simulation of a simplified postbounce stellar core. 

Differential rotation develops naturally in the outer core during stellar collapse 
[141]. Hence, such low-/3 instabilities could - in principle - occur in any PNS with (1) 
angular velocities at the edge of the inner core that are comparable to the pattern speed 
of unstable modes and (2) with sufficient total rotational energy and strong enough 
differential rotation to provide power for significant nonaxisymmetric deformation. 

The right part of figure 4 depicts the GW signal of a 20-Mo model of 
Ott et al. [115, 119] that experiences the low-/3 nonaxisymmetric instability. Shown 
are the + (top panel) and x (bottom panel) GW polarizations for observers located 
along the polar axis (red) and observers located in the equatorial plane at zero 
azimuth (black). The PNS and postshock region stay axisymmetric through bounce 
and exhibit GW emission only in the + polarization and away from the poles. The 
nonaxisymmetric dynamics reach sizable amplitudes around 30-40 ms and their m = 2 
components emit GWs in bar-like fashion with correlated strong emission in and 
hx that is strongest for observers located perpendicular to the equatorial plane. The 
simulation on which figure 4 is based was tracked to ~ 90 ms after core bounce, but 
the instability could potentially continue for hundreds of milliseconds until angular 
momentum has been redistributed by the instability or by another process (e.g., via 
the magneto-rotational instability [MRI]; see, e.g., [164, 165] and references therein). 

The simple and rough estimate of equation (1) for a slightly deformed PNS core 
rotating with a period of 2 ms matches the GW signal of the low-/3 instability shown 
in figure 4 quite well and other models of Ott et al. and Scheidegger et al. to within 
a factor of a few. Concretely, these studies found maximum GW signal amplitudes 
in the range of ^ 1-5 x 10^^^ (at 10 kpc) and total emitted GW energies Eqw of 
~ 5-15 X 10~* Mq c^. Most of the energy (> 50%) is emitted in a frequency interval 
of A/50 ~ 50-200 Hz about the frequency /peak at which dEQ^/df peaks, which is 
~ 900-940 Hz in the models of Ott et al. and Scheidegger et al. Comparing the just 
stated quantitative results with those summarized for rotating collapse and bounce in 
table 1, one notes that the GW emission due to nonaxisymmetric rotational dynamics 
does not necessarily dominate in signal amplitude, but does so clearly in total emitted 
energy. Its i?GW scales with the number n of "emission cycles," which, of course, 
depends on the duration of the nonaxisymmetric dynamics. 

As demonstrated by [96, 115, 119, 163], the potential for a strong enhancement 
of the GW signature of rotating core collapse by a low-/3-type rotational instability is 
great. However, the number of 3D core-collapse models in which such an instability 
was observed is still very limited and the systematic dependence on f3, degree of 
differential rotation, and thermodynamic and magneto-hydrodynamic configuration of 
the PNS remains to be established. Furthermore, it will be necessary to understand 
the long-term behavior of the instability and its interaction and competition with 
other shear instabilities such as the MRI. 
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6. Postbounce Convection and SASI 

Convective instability is a central feature of the postbounce evolution of core-collapse 
SN and has been discussed extensively (e.g., [2, 8, 21, 22, 89, 166-171] and references 
therein). According to the Schwarzschild-Ledoux criterion [172], convective overturn 
develops in the presence of negative radial entropy or composition gradients (see, e.g., 
[8, 167]). 

Convective overturn in postbounce SN cores is expected as entropy- and Icpton- 
gradient driven prompt convection which may occur immediately after bounce. Icpton- 
gradient driven PNS convection, and neutrino-driven convection in the postshock 
heating region. Convection will occur at postbounce times in virtually all core-collapse 
events, but can be weakened and limited to polar regions by positive specific angular 
momentum gradients in rapidly rotating cores [83, 89, 141, 173]. 

The standing-accretion-shock instability (SASI; see, e.g., [11, 21, 174-176] and 
references therein) is believed to be caused by either an advective-acoustic [11, 174] 
or a purely acoustic [175, 176] feedback cycle, leading to the growth of perturbations 
in the stalled shock. In axisymmetry, the dominant SASI modes are of £ = {1,2} 
character, while in the largely unexplored 3D case, power may go into azimuthal m 
modes. This can reduce the saturation amplitudes in the £ modes (see [177] for first 
3D results, but also [178] for a perturbativc analysis without symmetry constraints). 

The SASI grows to non-linear amplitudes over a timescalc of typically 200- 
300 ms after bounce with some dependence on initial accretion shock radius, neutrino 
luminosity, and accretion rate (e.g.. [10, 11]). Rapid rotation may delay the growth 
of the SASI and weaken the overall phenomenon in 2D [83] , while favoring azimuthal 
SASI modes in 3D [178]. 

Once active at significant amplitudes, the SASI heavily distorts the postshock 
region, affecting convection to great extent. In addition, it enhances and modulates 
accretion funnels that appear as kinks in the shock and channel low entropy material 
at high rates onto the PNS core [21]. 

Both convection and SASI are intrinsically multi-D phenomena and lead 
generically to time-varying mass quadrupole moments, hence, emit GWs. The 
emission of GWs by convection was first considered by Miiller and Janka [124] in 2D 
and 3D postbounce models with an approximate treatment of neutrino heating and 
cooling. Subsequent studies with more detailed neutrino transport and microphysics 
were carried out by Miiller et al. [97], Fryer et al. [179], Ott et al. [23, 42], and very 
recently by Marek et al. [98] . Kotake et al. [85] were the first to specifically address GW 
emission associated with the SASI using a simple neutrino heating/cooling scheme. 

In the following, we discuss the GW emission from convection and SASI, present 
new results that were performed with the VULCAN/2D SN code [22, 41, 180] and 
contrast them with previous work. VULCAN/2D is an axisymmetric Newtonian multi- 
group flux-limited diffusion (MGFLD) radiation-hydrodynamics code. It employs an 
unsplit arbitrary Lagrangian-Eulerian (ALE) scheme that allows for a central pseudo- 
Cartesian grid that transitions to a polar-type grid at a transition radius (see figure 4 
of [41]). 
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Figure 5. GW burst signal {h^D in units of cm) due to 2D prompt postbounce 
convection in a nonrotating 15-Mq model that employed the precollapse profile 
of [181]. The model calculations were performed with the VULCAN/2D code 
[22, 41, 180]. VULCAN/2D uses a central pseudo-Cartesian grid that is matched 
to a polar-type grid at a transition radius (20 km in sl5WW95HR and 30 km in 
sl5WW95). The transition introduces artificial perturbations that act as seeds for 
prompt convection. Model sl5WW95HR is a high-resolution variant of sl5WW95 
with ~ 45% more zones in the central region. The resolution - in other words - 
seed perturbation dependence of prompt convection is apparent. The GW signal 
data shown here are available from [155]. 

6.1. Prompt Convection 

As the stalling bounce shock passes through outer core material, it leaves behind 
a negative entropy gradient. Furthermore, following neutrino shock breakout and 
the associated burst of electron neutrinos (e.g., [8, 171, 182]), a negative lepton 
gradient arises at the outer edge of the PNS immediately below the neutrinosphere§. 
The two negative gradients lead to a convectively unstable region according to the 
Schwarzschild-Ledoux criterion [98, 183-185]. 

Neutrino losses and, to a limited degree, neutrino energy deposition behind the 
stalling shock, smooth out the large negative entropy gradient in the immediate 
postshock region, but prompt convection can still develop rapidly and last for > 10- 
20 ms if significant seed perturbations are present in the immediate postbounce flow 
[41, 179, 183-185]. 

The magnitude and distribution of seed perturbations in the central regions of 
iron cores in nature is presently unknown. However, it is not unlikely that the late 
burning stages of stellar evolution result in inhomogeneities in the iron core that do not 
smooth out completely between the end of core burning and core collapse (see also the 
discussion of large-scale precollapse asymmetries in section 9.2). Such inhomogeneities 
are frozen in during collapse and may act as seed perturbations for prompt convection 
[179, 186, 187]. 

Here, we present results from VULCAN/2D simulations using a nonrotating 
15-M0 progenitor model of [181]. Model sl5WW95 uses the standard simulation 
setup presented in [22, 23], while model sl5WW95HR is set up with 45% higher 

§ The neutrinosphere can be defined as the surface at which the optical depth for neutrinos , given 
by Ti, = f^dr\i/~^, is equal to 2/3. Here, Xu is the (energy-averaged) neutrino mean- free path. 
Note that the neutrino sphere is strongly dependent on neutrino energy and species. 
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radial and angular resolution inside ~ 300 km. The Cartesian-polar transition of 
the VULCAN/2D grid introduces artificial perturbations that converge away with 
increasing resolution and decreasing transition radius. These numerical perturbations 
act as seeds for prompt convection. 

Figure 5 depicts the GW signal of prompt convection as seen in the VULCAN/2D 
simulations. In both calculations, prompt convection is strongest inside a radius of 
60 km from the origin. However, one should keep in mind that the region in 
which prompt convection is most pronounced depends on the location and strength 
of the negative entropy gradient and on the location and magnitude of the seed 
perturbations. In the standard-resolution variant, the large seed perturbations lead 
to strong prompt convection and a correspondingly large GW burst setting in within 
~ 1 ms of core bounce. The maximum GW amplitude is ~ 6 x 10~^^ (at 10 kpc) 
and within ~ 20 ms, an energy of ^ 1.5 x 10"^" Mqc^ is emitted in GWs. The peak 
of the GW energy spectrum is at /peak ~ 680 Hz and 50% of the energy is emitted 
at frequencies 5/.5o = 150 Hz centered about /peak- The higher resolution calculation 
sl5WW95HR has considerably smaller seed perturbations, hence prompt convection 
is much weaker and and a smaller GW signal is emitted with |/irnax| ~ 9 x 10"^"^, 
£;gw » 6 X 10-12 Mqc2, /peak « 430 Hz, and A/50 » 250 Hz. 

The above numbers and the GW signals shown in figure 5 demonstrate the general 
seed-perturbation dependence of prompt convection and of the resulting GW burst 
signal. However, they constitute only two examples and do not bracket the parameter 
space of possible immediate-postbounce configurations. 

Marck et al. [98] also observed prompt convection in their very recent study 
focussing on the EOS dependence of neutrino and GW emission in the postbounce 
phase of a nonrotating 15-M0 model. They compared a model run with a soft variant 
of the Lattimer-Swesty EOS [140] with a counterpart model run with the rather stiff 
Wolff EOS [188]. Their models produce GW bursts associated with prompt convection 
that last ^ 30 — 40 ms and produce maximum amplitudes of ~ 1 — 2 x 10^^^ (at 10 kpc) 
with most of the emission occurring between ~ 60 and ~ 150 Hz. Furthermore, 
Marek et al. found that the model employing the stiffer nuclear EOS yields up to a 
factor of 2 larger maximum amplitudes and somewhat higher frequency emission than 
its softer-EOS counterpart. In the VULCAN/2D models discussed in the above, the 
GW emission from prompt postbounce convection occurs primarily at 400 — 700 Hz. 
This discrepancy with Marek et al. is most likely due to different numerical seed 
perturbations and different locations of the former and of the convectively unstable 
regions. 

For a more complete picture and an understanding of the possible prompt- 
convection GW signals, it will be necessary to conduct a comprehensive study of the 
quantitative systematics of prompt convection with variations in the magnitude and 
position of seed perturbations. Ideally such a study should be carried out with a code 
that uses a nearly perturbation-free grid and with clearly quantifiable perturbations 
added to the flow shortly before core bounce. 

6.2. PNS convection 

Owing to a negative radial lepton gradient, PNSs are unstable to convective overturn 
in a radial interval from ~ 10 km to ^ 30 km [8, 169, 171]. Convection sets in ^ 20- 
50 ms after bounce and may last for several seconds as the PNS slowly contracts and 
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Figure 6. GW signal due to PNS convection in the axisymmetric PNS model of 
Kcil ct al. [169, 189] as analyzed by Miiller et al. [97] in terms of the quadrupole 

pure-spin tensor harmonic amplitude A2Q, and h^-D = ^■\J~^ siri^ BA^q . The 
matter GW signal and the GW signal due to anisotropic neutrino emission are 
shown (for a discussion of the latter, see section 8). This figure is figure 5 of [97] 
and used by kind permission from the authors. 

deleptonizcs [169] after a successful SN explosion. If the explosion is not successful (or 
weak) and a black hole is formed, PNS convection and the associated GW emission 
stop abruptly. Also, if present, large-amplitude PNS core (;-modes can distort PNS 
convection and themselves lead to GW emission much stronger than that due to the 
convective motions (see section 7.1). 

Figure 6 displays the GW signal of PNS convection as analyzed by 
IVIiiller et al. [97] for the axisymmetric PNS convection simulations of Keil et al. [169, 
189] who followed the evolution of an isloated nonrotating PNS core from shortly 
after core bounce to 1.2 s. IVIiiller et al. reported typical strain amplitudes \h\ of 
2 — 5 X 10"^'^ (at 10 kpc), an energy emission Eqw of 1.6 x 10~^° Mqc^ over 1.2 s and 
a broad spectrum peaking around ~ 1300 Hz and with most of the energy emitted at 
700-1500 Hz. 

We analyze the GW signal emitted by PNS convection in the inner ^ 50 km 
of the aforementioned VULCAN/2D simulations sl5WW95 and sl5WW95HR in 
the interval of 20-250 ms after bounce. This interval is chosen, (1) to exclude 
as much as possible prompt postbounce convection, and, (2), because GWs from 
PNS core g-modes begin to dominate the convective GW signal at times later than 
~ 250 ms in the VULCAN/2D simulations (see section 7.1). For the two calculations, 
we find comparable GW strain amplitudes \h\ in the range of 1 — 4 x 10~^^ (at 
10 kpc). For the high- resolution variant specifically, we find an emitted energy i?GW 
of ~ 7.0 X 10~^^ Mqc^ with a peak of cLEq^ /df at ^ 350 Hz and most of the energy 
being emitted within ±150 Hz about this peak. By scaling i^cw to the 1.2 s considered 
by [97], we obtain ^ 3.4 x 10~^^ Mqc^ which is a factor of ~ 5 smaller than their 
value. The reason for this discrepancy is most likely the higher average frequency 
of the GW emission in the model considered by [97]. Note, however, based on the 
inset plot shown in figure 6 (figure 5 of [97]), one can estimate a dominant emission 
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frequency around ~ 300 Hz in the first ~ 100 ms of PNS convection. Thus, the 
higher-frequency emission must take place at later times. This is consistent with the 
fact that the PNS becomes increasingly compact with time, leading to convection at 
decreasing radii and generally on smaller spatial scales. As a consequence, one may 
expect a secular "chirp" of the GW emission frequency from ^ 300 Hz to > 1000 Hz 
over 1 s. 

6.3. Neutrino-driven Convection and SASI 

Neutrino heating below the stalled shock peaks in the inner part of the gain region:]: and 
decreases outward (see, e.g., figure 9 of [83]). This establishes a negative radial entropy 
gradient and makes the gain region unstable to convective overturn (e.g., [8, 124, 168] 
and references therein). Convection develops within 30-50 ms after bounce in the 
gain region which extends at these times typically from ^ 50-80 km out to almost the 
radius of the stalled shock at ~ 150-250 km. 

The left panel of figure 7 shows the GW signal of neutrino-driven convection and 
SASI in the VULCAN/2D postbounce calculations sl5WW95 and sl5WW95HR. The 
former calculation is carried out to ~ 850 ms after bounce, while the computationally 
more demanding higher-resolution latter is carried out to ~ 400 ms. In order to 
exclude PNS convection and PNS g-modes, only regions outside a spherical radius of 
60 km are taken into account in computing the GW emission. Since both convection 
and SASI occur in the same spatial domain and since the SASI-related distortions 
of the postshock region dynamically modify convection, the GW signals of these two 
hydrodynamic instabilities cannot be separated cleanly in a SN simulation. Also note 
that the GW emissions of PNS and neutrino-driven convection/SASI generally interact 
and in linear theory (and certainly in the Newtonian quadrupole approximation) 
superpose linearly. 

In the calculations considered here, the SASI reaches the non-linear regime around 
~ 300 ms and the shock and postshock flow exhibit strong temporally varying 
distortions at small and large scales (see, e.g., [21]). This is reflected by the GW 
signal in the form of more rapid time variation and a more pronounced modulation 
of the late-postbounce emission. For simulation sl5WW95, we find maximum GW 
signal amplitudes around 6 x 10~^^ (at 10 kpc; typical amplitudes are about a factor 
of 2 smaller), i?GW ~ 7.5 x 10"^^ Mq (for the entire 850 ms). The GW emission is 
rather broadband with most of the energy emitted between 100 and 500 Hz (though 
with a significant high-frequency tail which is emitted predominantly at postbounce 
times greater than ~ 400 ms). This is also reflected in the characteristic GW strain 
spectrum [191], 



shown in the right panel of figure 7 and contrasted there with initial and advanced 
LIGO rms noise levels. Note that /ichar of the higher-resolution calculation is globally 
lower due to the shorter period of postbounce time covered by this calculation. Such 
differences are not limited to simulations: After the onset of explosion, neutrino- 
driven convection is expected to cease quickly (within 50-100 ms [97]; note that PNS 

J The gain region is defined as the region of postbounce space in which neutrino heating dominates 
over neutrino cooUng. 
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Figure 7. Left: GW signal (h+ D in units of cm) due to postbouncc neutrino- 
driven convection in tlie core of a nonrotating 15- Mq presupernova model of [181] . 
Tile first 20 ms of postbounce evolution are cut in order to exclude the initial burst 
due to prompt convection. The standard-resolution calculation is run to ~ 850 ms 
after bounce. The development of the SASI to non-linear amplitudes around 
300 ms after bounce is reflected in increased amplitudes and more pronounced, 
higher-frequency variations. In order to study the resolution dependence of the 
convective GW signal, the high-resolution variant is run to ~ 400 ms after 
bounce. Because of the stochastic nature of turbulent convection, one cannot 
expect pointwise convergence of the GW signal, but rather rough agreement in 
signal amplitudes and characteristic frequencies. The latter is true for sl5WW95 
and sl5WW95HR. The waveforms are available for download from [155]. Right: 
''■char(/) spectra (equation 2) obtained from the sl5WW95 and sl5WW95HR 
waveforms for D = 10 kpc and compared to initial and burst-mode advanced 
LIGO /irms noise curves [190]. Note the broadband character of the convective 
GW emission and the high-frequency tail that is emitted predominantly at 
late postbounce times and, hence, is particularly pronounced in the standard- 
resolution variant that is run to ~ 850 ms after bounce. sl5WW95HR has overall 
smaller /i^har simply because it is run to only ~ 400 ms after bounce. 



convection continues). Hence, a real-life SN that explodes at earlier postbounce times 
will have a smaller integral and frequency-differential i?GW (and /ichar) from neutrino- 
driven convection/SASI than a counterpart that experiences a later or no explosion at 
all. The GW signals of sl5WW95 and sl5WW95HR agree well in average amplitudes 
and temporal variations in the first ^ 400 ms shown in the left panel of figure 7. For 
the long-term evolution of a turbulent system one generally cannot expect pointwise 
convergence for different resolutions. However, agreement "on average" of the GW 
signals is to be expected. 

Miiller et al. [97] analyzed the GW emission in 2D simulations of a nonrotating 
11.2-Mo and a slowly rotating 15-Mq progenitor. They followed the postbounce 
evolution of their models for 200-250 ms, included an approximate GR potential 
(see §3), Boltzmann neutrino transport in a ray-by-ray fashion [192], a relatively 
soft variant of the Lattimer-Swesty EOS [140] and comparable resolution to the 
aforementioned VULCAN/2D simulations. 

The average GW signal amplitudes and time variations found in the VULCAN/2D 
simulations agree roughly with those of IVIiiller et al. up to ^ 150 ms after bounce, 
but do not show the large enhancement (by a factor of up to ^ 20 and with rapid 
time variation around 500 — 800 Hz) of the GW emission that appear at later times in 
their simulations. IVIiiller et al. attributed this enhancement to the increased strength 
of convection, but there may also be a connection to the SASI excursions that become 
strong around these times in their models. Owing largely to that enhancement, 
Miiller et al. found that most of the GW energy is emitted in the interval of 500—800 Hz 
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and published total E'gw in the range of 10 — few x 10 ^ MqC^, roughly two orders 
of magnitude larger than seen in the VULCAN/2D simulations. 

The cause of the diserepancy between the Miiller et al. simulations and GW 
signals and those presented here is most likely related to differences in the employed 
EOS, the neutrino transport scheme and resulting neutrino luminosities and heating 
efficiencies, stalled shock radii, and in the development of neutrino-driven convection 
and the SASI. Furthermore, the enhancement of the GW signal may in part be due 
to PNS convection or PNS pulsations, since they did not exclude the entire PNS in 
their analysis. 

Kotake et al. [85] focussed on the SASI and employed a simplified approach for 
modeling postshock flow and neutrino transport. They assumed constant neutrino 
luminosities and starting from analytic postbounce conflgurations with a fixed 
accretion rate of 1 Mqs~^ that appear initially stable or only marginally unstable 
to convection, they studied the dependence of the SASI and SASI/convectivc GW 
emission on the assumed neutrino luminosity. The results of Kotake et al. confirm the 
expectation that higher neutrino luminosities lead to a more pronounced SASL more 
SASTinduced turbulence in the postshock region, and greater GW signal amplitudes 
and emitted energies. For the lowest luminosities considered {Li, — 55 Bs~^ each 
for I'e and 9e), they found maximum amplitudes of ^ 2.5 x 10"^^ and maximum 
hchar of ~ 10^^^ (dependent, of course, on the duration of the simulation) at /peak of 
100-200 Hz. The sum of the time-dependent ly^ and luminosities in the sl5WW95 
postbounce simulations with VULCAN/2D varies from 85 Bs^^ at ~ 100 ms to 
50 Bs^^ at > 400 ms after bounce. It is at all postbounce times smaller than the 
110 Bs~^ assumed by Kotake et al. Considering the dependence of the SASI and 
convective dynamics on the neutrino luminosity, it is not surprising that Kotake et al. 
found significantly stronger GW emission than the analysis of the VULCAN/2D 
simulations suggests. 

Marek et al. [98] extended the previously discussed work of Miiller et al. [97] and 
considered for a nonrotating 15-M0 model variations of the convective/SASI dynamics 
and the associated GW emission with the stiffness of the nuclear EOS. They followed 
their calculations to 400 ms after bounce and analyzed in detail the GW signal due to 
matter dynamics and anisotropic neutrino radiation fields (for a summary of the latter, 
see section 8.2). The qualitative features of the matter GW signals their models are 
qualitatively fairly similar to those obtained with VULCAN/2D (shown in figure 7). 
Marek et al. found that the GW signal amplitudes and typical emission frequencies 
increase with postbounce time and with the growing strength of the S ASI-driven shock 
distortions. Their models produce peak amplitudes around 5 x 10"^^ (at 10 kpc) and 
average amplitudes about a factor of 2 smaller. Furthermore, Marek et al. reported 
that the SASI-moderated asphcrical convective motions in the low-density postshock 
gain layer lead to GW emission at frequencies up to ~ 200 Hz. Those GWs, however, 
make up only the subdominant, low-frequency part of the total emission, which in 
their models is dominated by the GWs emitted from deeper, higher-density regions 
where aspherical accretion downstreams are decelerated and perturb the PNS core. 

Comparing a model run with a soft variant of the Lattimer-Swcsty EOS [140] with 
a counterpart model run with the stiff Wolff EOS [188], Marek et al. found that the 
Lattimer-Swesty EOS model yields higher-amplitude (by up to a factor of 2), higher- 
frequency (spectral peak at 600 — 800 Hz vs. 300 — 600 Hz) GW emission. These 
differences are a consequence of the more compact PNS and resulting harder neutrino 
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Table 2. Semi-quantitative summary of the GW emission by aspherical fluid 
motions associated with convection and SASI. The numbers for prompt convection 
are based on VULCAN/2D simulations presented here and in [23] and on 
simulations carried out by Marek et al. [98]. The GW emission characteristics 
of PNS convection and neutrino-driven convcction/SASI are based on the results 
of [23, 97, 98] and the VULCAN/2D simulations presented in this article. We 
provide estimates for the typical GW strain at 10 kpc, the typical emission 
frequency /, the duration of the emission At, and the emitted energy -Bgw in 
GWs. In addition, we list processes and factors that may limit the duration of 
the GW emission. All numbers, including various upper and lower limits, should 
be regarded as estimates that require confirmation by future studies and that, at 
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spectra and increased neutrino luminosities which, in turn, result in increased neutrino 
heating and more vigorous convection/SASI. This quantitative result is of particular 
importance, since it demonstrates that - despite the stochastic, untemplateable nature 
of the convective/SASI GW signal - key microphysical aspects can still be constrained 
via the observation of GWs emitted by convection/SASI. 

To conclude this section, we point out that in consideration of the still small set of 
mostly axisymmetric studies that have addressed GW emission from neutrino-driven 
convection/SASI, the currently available GW signal estimates should be regarded 
only as examples that may guide expectations. They are not yet robust and many 
more simulations will be needed to systematically (qualitatively and quantitatively) 
study the range of possible postbounce evolutions and corresponding incarnations of 
convective overturn. In particular, it will be important to quantify the dependence on 
progenitor structure and explore the dynamics of convection and SASI in 3D. 

Table 2 provides an overview of the GW emission from prompt, PNS, and 
neutrino-driven convection / SASI. We provide rough values for typical \h\, typical 
emission frequencies, emission durations, emitted energies, and list factors and 
processes that may inhibit the development or shorten the duration of the individual 
emission processes. 

7. Non-radial PNS Pulsations 

PNSs and cold NSs can pulsate in a multitude of ways (e.g., [61, 146]) and their 
non-radial pulsation modes of quadrupole and higher order emit GWs. Polar (even- 
parity) fluid modes are grouped into the lower-frequency (^-modes (buoyancy modes, 
restoring force gravity, require compositional or thermal gradients) and the higher- 
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frequency p- modes (restoring force pressure), separated in frequency by the /-mode 
(fundamental mode without nodes inside the star). Axial (odd-parity) fluid modes 
are modes whose restoring force is the Coriolis force, hence they require rotation 
and are degenerate at zero frequency in the nonrotating limit. Modes of mixed axial 
and polar character, so-called hybrid modes, exist as well (e.g., [146, 193]). r-modes 
(Rossby waves), a particular group of axial inertial modes, are strongly unstable to the 
GW-back-reaetion Chandrasekhar-Friedman-Schutz (CFS) instability [194, 195] at all 
rotation rates [61, 196]. In PNSs, r-modes have growth times of the order of seconds 
[197] and appear to saturate at low amplitudes [198]. Hence, they are unlikely to be 
relevant sources of GWs in the first second of a PNS's life. 

In addition to fluid modes there exist polar and axial modes of spacetime 
curvature, low-frequency s-modes and high-frequency w-modes, that are weakly 
coupled to matter (e.g., [60, 146] and references therein). 

Studies of non-radial pulsations of cold NSs as emission processes of GWs have a 
long history (see, e.g., [60, 61, 199, 200] and references therein). In the core-collapse 
SN context, non-radial PNS oscillation modes and damping by GW emission and 
other processes were first considered by MeDermott et al. [201, 202], Finn [203], and 
Reisenegger and Goldreich [204]. 

Recently, in a series of papers, Ferrari and collaborators [205-208] studied in GR 
perturbation theory the PNS non-radial mode spectrum and the GW emission from 
PNS pulsations on the basis of detailed spherically-symmetric PNS cooling models. 
They considered mode amplification by the CFS instability in rotating PNSs, but 
did not study the excitation of non-radial modes in the nonrotating limit. Assuming 
the presence of pulsational energy, they provided estimates in [205] for the typical 

I ft- 1 and i?GW necessary for a particular mode to be detectable by current and future 
GW observatories. According to their estimates for a nonrotating or very slowly 
rotating PNS, the most efficiently radiating mode, the lowest-order quadrupole g-mode 
(frequency / w 500 — 600 Hz), would need to emit GWs at a typical \h\ of ~ 8 x 10~^^ 
(at 10 kpc) and emit a total i?GW of ~ 4 x 10^® Mq to be detectable by initial- 
LIGO-class detectors at a distance of 10 kpc. In the ease of more rapid rotation, the 
(;-mode frequencies are lower and closer to the maximum detector sensitivity, leading 
to lower required \h\ and i?GW- 

Passamonti et al. [209, 210] explored the excitation of higher-order non-radial 
polar modes by non-linear couplings of radial and non-radial eigenmodes that are 
excited in the early postbounce ring-down phase of any core-collapse event§. Such 
couplings were first suggested by Eardley [47]. Passamonti et al. employed a GR 
perturbative framework and approximated the PNS by a simple polytropic stellar 
structure (which excludes (?-modes). Their results indicate that the radial and non- 
radial fundamental oscillation modes present in the very early postbounce phase will 
couple and give rise to a rich spectrum of daughter modes || . The strongest of these non- 
linear modes could be detectable by advanced LIGO-class detectors from a galactic 
event and their spectral information would provide additional insight into the high- 
density nuclear EOS and the internal structure of a newborn PNS. 

The pulsational mode structure of rapidly rotating NS/PNSs cannot be 

§ In a related work, Passamonti et al. [211] also studied the coupling of radial modes with non-radial 
axial modes. 

II Such radial-non-radial mode couplings have recently been observed also in numerical models of 
phase-transition-induced collapse of NSs [212]. 
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determined reliably by perturbative methods and GR hydrodynamic simulations 
must be carried out to obtain the spatial and spectral mode characteristics. 
Dimmclmeier et al. [199] recently carried out the first extensive study of the 
axisymmetric mode structure of rapidly rotating polytropic NS models in CFG 
GR [100], going beyond the Gowling approximation employed by most previous studies 
(e.g., [213] and references therein). The mode structure of rapidly rotating finite- 
temperature PNSs remains to be studied. 

7.1. Protoneutron Star g- Mode Pulsations in the Context of the Acoustic Mechanism 
for Core- Collapse Supernova Explosions 

Burrows et al. [21, 22], in their simulations that lead to the proposition of the acoustic 
mechanism for core-collapse SN explosions, have discovered the excitation of PNS 
core g-modes by turbulence and by accretion downstreams through the unstable and 
highly- deformed stalled shock that experiences the SASI. The pulsations damp by 
the emission of strong sound waves and do not ebb until accretion subsides. As 
discussed in [21, 22], these PNS core g-mode pulsations reach nonlinear amplitudes 
and act as transducers for the conversion of accretion gravitational energy into acoustic 
power that is deposited in the postshock region and may be sufficient to drive an 
explosion. Although promising, this acoustic SN mechanism remains controversial 
(see, e.g., [9, 24, 25]) and has yet to be confirmed by other groups. In particular, 
Weinberg et al. [25] pointed out that the saturation amplitude of the PNS g-modcs 
may be limited by non-linear model coupling via a parametric instability that cannot 
be resolved in present simulations. In addition, there are indications [17, 83] that 
rapid rotation significantly weakens convection and modifies the SASI and thus may 
weaken the driving mechanism for the PNS core (7-modes. 

Burrows et al. found that the fundamental £ = 1 g-mode of the PNS core is most 
easily and first excited in their models. Higher-order eigenmodes and, through non- 
linear couplings, harmonics of eigenmodes and other modes with complicated spatial 
structures appear at later times. Ott et al. [23] analyzed the GW emission from the 
quadrupole components of the PNS core pulsations in the original nonrotating 11-Mq 
model presented in [21] and, in addition, for a slowly-rotating 15-M0 model and a 
nonrotating rather extreme 25-71/© model with a very extended and massive iron core. 
For these models they found maximum GW amplitudes |/imax| of ^ 1.3 x 10"^"'^ (at 
10 kpc) for their ll-M© and 15-M0 models and - 5 x 10"^" (at 10 kpc) for the 
25-M0 model. The total emitted -EgwS were around 1.5 x 10"^ Mq c? in the ll-M© 
and 15-M0 model, while the 25-M0 model emitted an amazing 8 x 10~^ Mq (? . GW 
back-reaction effects, though neglected in [23], are probably relevant for this latter 
model [42]. Most of the emission in these models took place in frequency space at 
~ 650 ± 100 Hz and - 900 ± 100 Hz. 

In [22], Burrows et al. considered an extended set of nonrotating progenitors, 
covering the mass range from 11.2 M© to 25.0 M© solar masses, and using 
presupernova models of [137, 214]|. In the following, we present their previously 
unpublished GW signals and summarize the key signal characteristics in table 3. The 
waveforms and E'gw spectra of all models can be obtained from [155]. 

The left panel of figure 8 depicts the GW signal obtained from models sl5.0 

X Note that the precollapse models used in [21, 23] were taken from [181] and [142]. Model s25.0 of 
[22] should not be confused with model s25WW of [23] . 




Figure 8. GW emission in a set of representative models from [22]. Left: 
GW wave signal (fe+D in units of cm) as a function of postbounce time in the 
nonrotating models s25.0 (red) and sl5.0 (green). The burst shortly after core 
bounce is due to prompt convection (see section 6). Around 300-400 ms, the 
GW emission becomes dominated by the growing g-mode amplitudes, but stays 
at relatively small amplitudes initially. The core pulsations become more vigorous 
at later times, leading to the emission of large-amplitude GWs. Explosion sots in 
around ~ 1.0-1.4 s after bounce and in some cases is unipolar (sec [21, 22]) - the 
star explodes on one side while accretion continues on the other. All waveforms 
are available from [155]. Right: Characteristic GW strain spectra (equation 2 
and [191]) of various models of [22] contrasted with the initial and advanced LIGO 
rms noise curves [190] for a source distance of 10 kpc. Note that the spectra of 
all models exhibit a narrow primary peak near 950 Hz and a broader secondary 
peak around 600-700 Hz. 



(15 Mq) and s25.0 (25 Mq) as a function of postbounce time. There is little variation 
in the qualitative features of the GW signal throughout the model set considered 
by [22]. Hence, models sl5.0 and s25.0 are quite representative picks. Core bounce 
in these models is marked by a GW burst that is due to prompt convective overturn 
setting in very shortly after bounce (see section 6.1). Typical peak GW amplitudes are 
around ^1x10"^^ at 10 kpc and most of the energy is emitted in a frequency interval 
of 500-800 Hz. The prompt convection subsides quickly, and during the interval from 
~ 20 ms to ^ 300 ms the GW signal is dominated by convection in the PNS and in 
the postshock region. The PNS core pulsations reach significant amplitudes at 300- 
400 ms after bounce and the GW emission from their quadrupole components starts 
to dominate the overall GW signal. From ~ 300 ms to ~ 1.1 s the models emit quasi- 
continuous GWs with maximum amplitudes in the range of ~ 3 x IQ-^^-S x lO"^'^ 
at 10 kpc, modulated on a timescale of ^ 100 ms by variations in the SASI and 
accretion-dependent excitation (and damping) of the core pulsations. At postbounce 
times greater than ^ 1.1 s, the quadrupole components of the core pulsations reach 
very large amplitudes (possibly through the excitation of an ^ = 2 eigenmode [23]), 
leading to a significant increase in the GW amplitudes up to ~ 7 x 10"^^ (~ 4 x 10"^^) 
at 10 kpc in model s25.0 (sl5.0). The other models studied by Burrows et al. [22] reach 
comparable |ft-max| (see table 3). 

All models of Burrows et al. [22] explode between ^ 1.0 s and 1.4 s after bounce 
and the simulations were stopped when the shock reached the outer boundary of the 
computational grid. In models that exploded globally, the PNS core pulsations and 
the corresponding GW signal generally subsided after the onset of explosion. However, 
in a few models, the explosion was unipolar, setting in along one of the poles while 
accretion continued on the opposite side, sustaining to some extent the core pulsations 
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Table 3. GW data summary for models calculated by Burrows et al. [22] 
whose GW signals have not previously been published. In addition, the GW 
emission characteristics of the models of [21, 23] are listed for completeness. 
At is the amount of postbounce time covered in each model, |/i+,max| is the 
overall maximum GW amplitude scaled to 10 kpc distance, /ichar.max is 
the global maximum of the characteristic strain spectrum [191], /peak is 
the frequency of the global maximum and Sf is the FWHM. Data for the 
secondary pronounced peak in hchnvif) are given and all related quantities 
have the subscript 2. -Egw is the total energy radiated in GWs. Note that 
the values given for | ftmax j in [23] are slightly different from those presented 
here for models of [23] . This is due to a small error in the analysis routines 
used by [23]. Also note that model s25WW of [23] uses a rather extreme 
25-Mq precollapse model of [181] with a very extended and massive iron 
core and a very shallow density gradient. 
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and prolonging the GW emission. 

In the right panel of figure 8, characteristic strain spectra ft.char(/) (equation 2; 
for a source distance of 10 kpc) are shown for models sl3.0, slS.O, and s25.0. More 
massive progenitors generally tend to have higher accretion rates§, tend to explode 
later, reach higher | /imax | , have higher total E'gw , and more narrow peaks in diScw / df 
(see table 3). As a consequence, the maximum values of /ichar generally increase with 
progenitor mass. 

It is interesting to note that /ichar exhibits two pronounced peaks in all considered 
models. The primary, quite narrow peak is centered in frequency around 900-1000 Hz 
while the secondary, smaller and broader peak is located in the frequency range 500- 
700 Hz. This is due to the complicated time-dependent mode structure present in the 
pulsating PNS core and is also highlighted by figure 9 which provides a frequency-time 
analysis of the GW energy spectrum dE'cw/^/ of model sl5.0. Due to changes in the 
hydrodynamic and thermodynamic structure of the PNS, the g-mode frequencies vary 
(as do the frequencies of other mode families) with time [205]. In addition, higher- 
order harmonics and modes with complicated spatial structure are excited and their 
quadrupole parts emit GWs. At early to intermediate postbounce times, the emission 

§ Note that the presupcrnova models of [137] behave non-monotonically in the scaling of (iron) core 
mass and density profile in the mass range from ~ 15 to ~ 20 Mq. See also [22, 108]. 
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Figure 9. Frequency-time evolution of the GW energy spectrum dEQ-^/df in 
model sl5.0 of [22], computed with a 50 ms sampling interval and a step width 
of 1 ms. Shown is the interval from 0.4 s after core bounce to the end of the 
simulation. Given the complex thermodynamic structure of a hot PNS, the mode 
spectrum reflected in dEQtjv/df is rather complicated and exhibits significant 
temporal variations in mode frequencies and preferred modes. Note that the two 
strong emission regions around 900-1000 Hz and 600-800 Hz correspond to the 
the primary and secondary peak in /ichar(/) shown in the right panel of figure 8. 



occurs mainly from quadrupole components with frequencies around ^ 600 Hz and 
850 Hz. The frequencies initially show a weak upward trend, but decrease beginning 
at around ^ 700 ms. At late times, most of the energy is being emitted around 
~ 850 Hz and ^ 600 Hz-700 Hz, but this time by quadrupole components of modes 
who have descended from higher frequencies. 

The PNS core g-mode pulsations may arguably be so far the strongest proposed 
emission process for GWs in core-collapse SNe. Observing GWs due to PNS g- 
modes and capturing the time-evolution of the GW spectrum would provide us with 
invaluable information that could be used, via comparison with model calculations, to 
reconstruct the hydrodynamic and thermodynamic structure and evolution of a very 
young PNS. In addition, and since the acoustic mechanism for SN explosions and the 
strong PNS core (7-modes are almost invariably linked, the observation or the non- 
observation of the g-mode GW signal from a nearby core-collapse SN may be used to 
test the acoustic SN mechanism (see also section 11). 

Nevertheless, it is important to keep in mind that the viability of the acoustic 
SN mechanism is still unclear and that the dependence on dimensionality (2D vs. 
3D), rotation, equation of state, description of gravity (Newtonian vs. GR) and the 
possibility of mode saturation at low amplitudes [25] remain to be explored in detail. 

8. Anisotropic Neutrino Emission 

Not only asphcrical fluid motion, but in more general terms, any accelerated transport 
of energy with a non-zero quadrupole and/or higher-order component emits GWs. 
For the case of anisotropic radiation of neutrinos from a distant point source, this 
has been first realized in linear theory via the direct solution of the inhomogeneous 
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wave equation by Epstein [81] (but. see also Turner's independent derivation in the 
zero-frequency limit [82]). 

Burrows & Hayes [123] and Miiller & Janka [124] were the first authors to 
implement the formalism. It has since been employed in a number of other studies 
[23, 35, 85, 97, 98, 179, 215]. In the following, we present a short overview on the 
formalism as used by [123, 124] in axisymmetry. More details and the generalization 
to 3D can be found in [35, 124]. 

For axial symmetry, [123, 124] write the dimensionless GW strain for an observer 
positioned in the equatorial plane as 

hll^{t) = — / a{t')L,{t')dt' , (3) 



where L^(t) is the total neutrino luminosity and a(t) is the instantaneous neutrino 
radiation anisotropy that includes the transverse-traceless projections [124]. It is 
defined as 

where dL,j{VL^t) / d£L is the energy radiated at time t per miit of time and per unit of 
solid angle into direction with 

\E'(i9, (ys) represents the angle dependent factors in terms of source coordinate system 
angles i? and <p and depends on the particular GW polarization and the observer 
position relative to the source. In axisymmetry, ~ everywhere, and = 
along the axis of symmetry. For an observer located in the equatorial plane, observing 
the + GW polarization, '^{■d.ip) is given [35, 124] by 

^ , „ , , . „ , cos^ d — sir? d sir? in 

vl/(z?,^) = (l + sm^cos^) . 2 c ■ 2 ■ (6) 

cos^ v + sm v sm (p 

In axisymmetry, there is no dependence of the luminosity. By integrating ^'(i?, (p) 
over equations 4 and 6 combine to [35, 85], 

Note that the GW signal due to neutrinos observed at time t + D/c contains 
contributions from anisotropics in the neutrino radiation field at all times prior t. This 
leads to a memory effect in the GW signal, leaving behind a constant ("DC") offset 
after the anisotropic neutrino emission subsides. Largely-aspherical mass ejection can 
lead to a similar GW memory (see section 9). The implications and detectability of 
such GW bursts with memory were discussed in [34, 216]. 

In the context of massive star collapse and core-collapse SNe, anisotropic neutrino 
radiation and the associated emission of GWs may arise (a) from rotationally-deformed 
PNSs [35, 42, 97, 215], (b) from convective overturn and SASI [23, 42, 85, 97, 98, 124], 
and (c) from global asymmetries in the (precollapse) matter distribution [123, 179]. 

The extraction of GWs due to neutrino radiation field anisotropics requires that 
the underlying simulation was carried out with some form of neutrino transport. 
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Unfortunately, the set of presently employed approaches to neutrino transport 
generally yield different results for the degree of anisotropy of the neutrino radiation 
field for the same hydrodynamic configuration. For example, any ray-by-ray scheme 
(Boltzmann transport or MGFLD) that implements 2D neutrino transport via many 
solutions of ID transport problems along radial rays, even if neighboring rays are 
coupled, tends to overestimate local and global anisotropics. On the other hand, a 
full-2D MGFLD approach tends to smooth out anisotropics in the radiation field, in 
particular at low optical depths. Ideally, neutrino transport should be carried out 
in 2D (or 3D) with energy and momentum-space angle-dependence. Ott et al. [83] 
have recently carried out such simulations using the S'„ method (e.g., [217]), but due to 
limited number of momentum-space angles that could be used in long-term evolutions, 
anisotropics in the neutrino radiation field are dominated by poor numerical resolution. 

When interpreting any results on the GW signal emitted by anisotropic neutrino 
radiation obtained with current codes, the dependence on the neutrino transport 
formulation and numerics should be kept in mind. 

8.1. Rapid Rotation 

Because of the action of centrifugal forces, a rapidly-rotating PNS has an oblate shape. 
The neutrino radiation field and its energy-dependent neutrinospheres follow the 
matter distribution. Because of the more compact polar and more extended equatorial 
density distribution, neutrinos decouple from matter at smaller radii along the poles 
than on the equator. Depending on the degree of rotational deformation, this can 
lead to very large pole-equator asymmetries in the neutrino radiation fields, resulting 
in the emission of GWs (see, e.g., [83, 215, 218, 219]). Since the qualitative shape of 
the asymmetry does not vary with time (viz. the PNS stays oblate), the GW signal is 
monotonically growing in amplitude and exhibits slow time variation. 

Dcssart et al. [215] and Ott [42] analyzed the GW emission due to aspherical 
mass motions and asymmetric neutrino radiation fields in the axisymmetric Newtonian 
accretion-induced collapse (AIC) simulations of Dcssart et al. [215]. An AIC event in a 
massive ONeMg white dwarf is expected to leave behind a PNS that is similar in many 
ways to a PNS formed in standard iron core collapse. The two AIC models considered 
by Dcssart et al. [215] yielded rapidly rotating PNSs with postbounce /?s of ~ 0.05 
and ^ 0.25 for their 1.46-Mq model and 1.92-Mo model, respectively. Dessart et al. 
employed a 2D MGFLD approach to neutrino transport and, hence, the numbers given 
below for the GW emission due to neutrinos may be underestimating the true signal 
strength. 

Figure 10 shows the contribution of anisotropic neutrino emission to the GW 
signal in the AIC models. At bounce, the signal shows the greatest variation due to the 
strong electron-neutrino burst. After bounce, the amplitude slowly but continuously 
rises. At the end of the simulations, [/imaxj at 10 kpc is ^ 3.6 x 10^^^ 0.55 x 10^^^) 
in the 1.92-M0 (I.46-M0) model. The amplitude continues to grow in the subsequent 
cooling phase, but its growth rate decreases continuously, since the neutrino luminosity 
is dropping at a greater rate than the anisotropy parameter a is increasing. In the 
frequency domain, the GW emission due to neutrinos occurs at low frequencies and 
dEQyj/df peaks in the range of 0-10 Hz. The total emitted i^cw is ~ 2 x 10~^^ Mqc^ 
(- 1 X 10-12 AIqc^) in the 1.92-M0 (I.46-M0) modcljl. 

II Note that the numbers for -Egw given here differ from those stated in [215]. This is due to an 
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Figure 10. GW signal from the rapidly-rotating accretion-induced collapse (AIC) 
simulations of Dcssart ot al. [215] as analyzed in [42]. The 1.92-A?fQ model has a 
postbouncc 13 of ^ 0.25 while the 1.46--Mq model is rotating much more slowly 
with postbouncc /3 < 0.05. Top: GW signal {h+D in units of cm) emitted by 
anisotropic neutrino radiation. The time is given with respect to the time of 
bounce and the first 500 ms after bounce are shown. Note the very slow time 
variation of the signal which is due to the time integral in equation 3 and to the 
slowly varying anisotropy of the radiation field owing to the secular contraction of 
the PNS. The GW signal due to neutrinos was extracted at a spherical radius of 
300 km (200 km) in the 1.92-M0 (I.46-M0) model. Bottom: Combined matter 
and neutrino GW signal. The GW signal data are available for download from 
[155]. 



Kotake et al. [35] showed example results of their unpublished work on the GW 
signal from rapidly-rotating iron core collapse. Unfortunately, they provided only the 
GW signal from core bounce and did not provide overall signal characteristics that 
could be compared with the AIC models of Dessart et al. 

Miiller et al. [97] studied the GW emission due to neutrinos in a relatively slowly 
rotating IS-Af© model. They found an almost monotonically-growing GW signal 
qualitatively similar to those shown in figure 10 with a maximum amplitude at 10 kpc 
of ^ 1.3 X 10"^^ at the end of their simulation 270 ms after bounce). Since their 
15-M0 model, in contrast to the AIC models of Dessart et al., exhibited significant 
convective overturn, its GW signal contains higher-frequency (~ 20 — 100 Hz) 
components that are due to either locally-enhanced neutrino emission from rapid 
downfiows or variations in neutrino absorption due to convective overturn. 

Given the above listed results, the GW signal due to anisotropic neutrino emission 
associated with rapid PNS rotation is emitted at too low frequencies and is not 
sufficiently energetic to be a good candidate for detection by current and planned 
ground-based GW observatories even if the stellar collapse event occurs within the 
Milky Way. The situation may be different in the case of moderate rotation studied 
by Miiller et al. [97] in which the GW signal contains higher-frequency contributions 

error in the GW analysis of [215]. 
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from convcction/down-strcam-induccd neutrino radiation anisotropics. According to 
the results of Miiller et al. , the neutrino GW signal of such cores may be detectable 
with advanced LIGO-class detectors for a galactic source. 

8.2. Convection and SASI 

Convection and SASI introduce asphericity into the postshock flow, leading, among 
other things, (1) to spatially varying neutrino cooling/heating in the region between 
shock and PNS core and (2) to accretion funnels that allow for rapid downflow 
of material and locally enhanced neutrino emission where an accretion funnel hits 
the PNS core. Janka & Miiller [124] first investigated (though with simplified 
neutrino transport and an artificial inner core boundary) the GW emission from 
anisotropic neutrino radiation fields owing to postbounce convection and accretion 
downstreams. Since they provided only GW signals for which neutrino and baryonic 
matter contributions were summed up, it is difficult to extract the neutrino component. 
A rough estimate gives modest peak amplitudes in the range of 5 x 10~^* to < 1 x 10^^^ 
at 10 kpc and the typical slow waveform variation. 

More recently, Miiller et al. [97], using a ray- by-ray Boltzmann neutrino-transport 
scheme and a gravitational potential whose monopole term is general relativistic, 
investigated the emission of GWs due to neutrinos in the first 200 — 270 ms after 
bounce in a slowly-rotating 15-Mq model (discussed above in section 8.1) and a non- 
rotating II-M0 model. In the latter, the GW signal due anisotropic neutrino emission 
reflects the stochastic nature of postbounce convection and varies slowly (typical 
frequencies around 10 — 20 Hz) with no clear long-term trend, reaching maximum 
amplitudes of ~ 1.4 x 10^^^ at 10 kpc. Due to the low-frequency variation, the 
emitted Eqw can be expected to be small, was not stated explicitly by the authors, 
and is likely to be in or below the ball park of the numbers mentioned in section 8.1 
for rapidly-rotating AIC models. In contrast to their slowly-rotating 15-Mq model, 
the neutrino GW signal of the non-rotating U-Mq model is unlikely to be detectable 
even by advanced LIGO-class observatories and when occurring at 10 kpc distance. 

In addition to the 15-Mq and 11-Mq models, Miiller et al. also considered 
the GWs emitted by anisotropic neutrino radiation in a PNS convection model of 
Keil et al. [169, 189] (sec section 6.2 and figure 6). The GW signal is qualitatively and 
quantitatively similar to what Miiller et al. observed in their full 11-A/q progenitor. 
However, before explosion, the true asymptotic neutrino anisotropy must be extracted 
outside the optically-semi-transparent postshock layer (see, e.g., [98]) which was not 
included in the PNS model of Keil et al. At post-explosion times, the neutrinos 
practically free-stream from the PNS surface and the GW signal from the PNS 
simulation becomes more relevant. 

Using 2D-MGFLD, Ott et al. [23, 42, 155] found for three long-term postbounce 
models in the mass range from 11 to 25 Mq maximum amplitudes [/imaxl of ^ 
1.3 X 10~^^ to ~ 5.5 X 10~^^ (at 10 kpc) built up over the ^ 1 s of postbounce 
time covered by the simulations and extracted at a radius of 200 km. The neutrino 
GW signals in these models exhibit a systematic trend to negative, in absolute value 
continuously growing amplitudes. The total emitted -Eqw due to neutrinos were 
small, below ~ 1 x lO^^'^ MqC^ and most of the energy emission took place at 
frequencies below ~ 10 Hz. Since MGFLD has the tendency to smooth out radiation- 
field anisotropics [83], the above numbers unfortunately depend on the extraction 
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radius and should be regarded as underestimates. 

Kotake et al. [85], using a simplified approach to neutrino heating and cooling, 
studied the GW emission from angular variations in the neutrino cooling in their 
models that were focussed on the non-linear development of the SASI. They found 
quasi- monotonically growing GW signals that showed slow variations (1-50 Hz) and 
reached maximum values after ^ 500 ms of 2.0-3.5 x 10~^° (at 10 kpc). The GW 
signal amplitudes found by Kotake et al. are systematically positive which the authors 
attributed to SASTenhanced polar neutrino emission. The emitted i^cw were in the 
range of 2 — 6 x 10~^° Mq(? . These large numbers are a most likely a consequence of 
the neutrino luminosities they imposed which were a factor of 5 to 10 larger than in 
more realistic simulations [23, 97] and lead to very vigorous overturn and SASI. 

Marek et al. [98] studied the GW emission from anisotropic neutrino radiation 
fields in a nonrotating 15-M0 model that was run with two nuclear EOS of 
different stiffness (see also the discussion of their study in section 6.3). Similarly to 
Ott et al. [23, 42], but in disagreement with the less detailed study of Kotake et al. [85], 
they found a continuously (in absolute value) growing negative GW signal. According 
to their detailed analysis, this negative signal is a consequence of enhanced equatorial 
neutrino emission in the postshock region. At 400 ms after bounce, their Lattimer- 
Swesty EOS [140] model reaches \h\ of ^ 1.2 x lO^^^ (at 10 kpc) while the model 
run with the stiffer Wolff EOS [188] exhibits a smaller \h\ of - 5.7 x 10"22_ 
\h\ are larger than the typical amplitudes of GWs emitted by the matter dynamics 
associated with convection/SASI (section 6.3). However, in agreement with previous 
studies, Marek et al. [98] provided GW spectra indicating emission predominantly at 
frequencies < 20 Hz, making the component of the overall GW signal associated with 
anisotropic neutrino emission very difficult to detect for ground-based detectors whose 
sensitivity is limited by gravity-gradient and seismic noise at low frequencies (e.g., 
[33]). 

8.3. Global Asymmetries 

Large-scale density perturbations that may be due to inhomogeneities in the iron core, 
silicon, and/or oxygen shells lead to globally asymmetric postbouncc mass motions 
resulting in large angular variations in the neutrino luminosity (see also section 9.2). 
The GWs emission due to neutrinos in such a globally asymmetric core-collapse event 
has first been considered by Burrows and Hayes [123] in 2D, using a ray-by-ray gray 
flux-limited diffusion approach [167]. Their single model with a 15% precollapse 
density perturbation near the north pole produced a GW signal with a large [/imaxj 
of ~ 3.5 X 10^^^ (at 10 kpc). The signal, shown in the right-hand panel of figure 11, 
exhibits a strong burst at core bounce (with ~ 500 Hz characteristic frequency) and a 
subsequent slow growth to the final amplitude (the model was tracked to only ~ 70 ms 
after bounce). The total emitted -Eqw is around '--^ 3 x 10"^" Mqc^ , the overall largest 
£'gw from anisotropic neutrino emission found for any presently published model. 

Fryer et al. [179], also using gray flux-limited diffusion, but employing a 3D 
Newtonian smooth-particle hydrodynamics (SPH) scheme, studied four models with 
density perturbations of 25% to 40% in 30-degree wedges about the north pole. They 
found GW signals due to anisotropic neutrino emission with maximum amplitudes 
of the same order of magnitude as those reported by [123]. However, their signals 
exhibited much less time variation and a weaker burst associated with core bounce. 
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Figure 11. Left: Grcy-scalc map showing tiio entropy distribution ~ 50 ms into 
the explosion of the 2D asymmetric collapse model of Burrows and Hayes [123]. 
The physical scale is 4000 x 4000 km and velocity vectors are superposed. This 
plot corresponds to figure 1 of [123] and is used by kind permission from the 
authors. Right: GW signal due to anisotropic neutrino emission (dashed line), 
matter (dotted line) and combined matter+neutrinos (solid line) in the same 
model of [123]. This plot corresponds to figure 3 of [123] and is used by kind 
permission from the authors. 

suggesting much lower energy emission (Fryer et al. did not provide values for i?Gw)- 

9. Other GW Emission Mechanisms 
9.1. Aspherical Outflows 

Core-collapse SN explosions are unlikely to be perfectly spherically symmetric and 
all currently considered explosion mechanisms rely on the breaking of symmetry 
IVIagnetically-driven explosions exhibit jet-like bipolar outflows (e.g., [17] and 
references therein), while explosions driven by the neutrino mechanism (e.g., [9]) 
or the acoustic mechanisms (e.g., [22]) can range from predominantly unipolar or 
bipolar outflows to nearly spherically symmetric explosions. In addition to the various 
explosion scenarios, precollapse large-scale asymmetries in the central regions of the 
star may also favor a largely asymmetric explosion. 

The GW emission from rapid aspherical outflows has been considered by 
Obergaulinger et al. [90, 125] in the context of 2D magneto-rotational core collapse 
in Newtonian gravity and approximate GR (see section 3) and by Shibata et al. [117] 
in 2D GRIVIHD. Both groups of authors employed polytropic precollapse models, a 
simple polytropic + ideal gas ("hybrid") EOS, and neglected neutrinos. At this level 
of approximation, the shock generally does not stall and a prompt explosion occurs 
after bounce. 

Obergaulinger et al. and Shibata et al. both reported a quasi-monotonically 
growing contribution to the GW signal from bipolar outflows, leaving behind a GW 
memory akin that discussed in the context of GW emission from anisotropic neutrino 
radiation flelds (see section 8). An example GW signal taken from the study of 
Obergaulinger et al. is shown in figure 12. Shibata et al. [117] pointed out that the 
GW signal systematics can be understood by considering an outflow of mass m in the 
z direction with slowly changing velocity. In this case, the contribution to the GW 
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Figure 12. GW signal of the Obergaulinger et al. [90] model A1B3G3-D3M13 
in terms of the quadrupole pure-spin Tensor harmonic amplitude j4^^, and 

h^-D = ^ sin'^ 6A2Q. The dashed line gives the contributions of regions 
with radius r < 60 km while the solid line shows the total signal. The postbounce 
offset due to rapid aspherical outflows in the outer regions is clearly visible in the 
total signal. This figure corresponds to the lower-left panel of figure 8 in [90] and 
is used by kind permission from the authors. 



signal by the outflow is |ft-out| oc 2mv^. Once the outflow has reached a quasi-steady 
state, the mass ejection rate can be assumed to stay roughly constant, hence the mass 
in the outflow increases roughly linearly [17, 117]. Using this, the appropriate factors 
of G and c, and imposing a scaling on m and Vz, one finds [117] 

i'--i^"™(o:S;)((n:)'-=- 

Similar to the case of anisotropic neutrino emission, the slow variation of |/iout,-i-| leads 
to emission of only small Eqw that are emitted predominantly at low frequencies. 

The approximate |/iout.+ | given by equation 8 is in rough agreement (assuming 
m — 0.1 Mq and Vz = 0.1c) with the GW memory amplitude of ^ 6.5 x 10^^^ (at 
10 kpe) that Shibata et al. found after ~ 15 ms of postbounce evolution in the most 
extreme model of their limited model set. Obergaulinger et al. [90, 125] performed a 
more extensive set of calculations, but did not publish the final GW signal memory 
amplitudes for all their models. Those that were shown in [90, 125] are in rough 
agreement with equation 8 and the results of Shibata et al. [117]. 

Finally, it is important to point out that the IVEHD-driven bipolar outflows 
observed in the microphysically-simple calculations of [90, 117, 125] occur very early 
(around ~ 5-20 ms) in the postbounce evolutions of their models. IVIore realistic 
models computed by Burrows et al. [17] (who did not analyze the GW emission in 
their models) require of the order of ^ 100 ms to launch the bipolar outflow. This 
longer delay translates to lower accretion and mass expulsion rates and and may result 
in smaller GW signal memory amplitudes than observed by [90, 117, 125] for their 
very early explosions. 
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9.2. Global Precollapse Asymmetries 

Prior to core collapse and after silicon core burning has ceased, nuclear burning 
continues in convectively unstable silicon and oxygen burning shells at the outer edge 
of the iron core (see, e.g., [137]). The large- and small-scale density, composition and 
thermal perturbations induced by convection in these layers and at the outer edge of 
the iron core could become large [187, 220-222], perhaps up to ^ 10%, and would 
essentially be frozen in during collapse [186]. 

The effect of these perturbations on the postbouncc hydrodynamics depends 
on their magnitude and distribution. If they occur primarily on small scales and 
reach deep into the iron core, they may act as seeds for prompt and neutrino-driven 
convection (see section 6 for the associated GW signals). Larger-scale perturbations, 
on the other hand, can lead to strongly aspherical explosions into the direction of 
lowest density. Such largely aspherical ejection of matter ~ as just discussed in the 
previous section 9.1 - can lead to the burst emission of GWs with GW memory (e.g., 
[216, 223]) and, in the case of primarily unipolar outflow, may be responsible for pulsar 
birth kicks (e.g., [224, 225] and references therein). 

Burrows & Hayes [123] were the first to perform a numerical study of core 
collapse with a globally asymmetric presupernova configuration. They employed an 
axisymmetric Newtonian code with gray flux-limited neutrino diffusion along rays. 
Their single model was a IS-M© presupernova configuration perturbed by a 15% 
reduction of the density of the outer core within 20 degrees of its north pole. The left 
panel of figure 11 displays the aspherical explosion obtained in their model. It erupted 
preferentially through the lower-density polar regions and lead to a recoil velocity of 
the PNS of ~ 530 km/s. The aspherical matter dynamics resulted in a GW burst signal 
(shown in figure 11) that set in ~ 5 ms after bounce and reached maximum amplitudes 
of ~ 1 X 10~^^ (at 10 kpc) at early times with rapid time variation around ~ 500 Hz 
and an asymptotic memory of ^ 1.5 x 10^^^ (at 10 kpc) built up over ~ 70 ms. The 
combined neutrino and matter GW signal exhibited a memory of ^ 2.3 x 10"^^ (at 
10 kpc) and the total emitted energy in GWs was ^ 1.1 x 10~^ AIq (? . 

Fryer et al. [179, 226], using their 3D SPH scheme with gray flux- limited diffusion, 
performed a set of four calculations of a 15-Mq presupernova star. Three calculations 
were perturbed with an initial 30% to 40% density reduction in the oxygen and silicon 
layers in a 30-degree wedge centered about one of the poles. One calculation was set 
up with 25% lower density in a 30-degree polar wedge throughout iron core, silicon and 
oxygen shells. Despite the larger perturbations present in their models. Fryer et al. did 
not find strongly asymmetric explosions and reported neutron star kick velocities below 
~ 200 km/s, arguing that momentum transfer by neutrinos was partly counteracting 
the hydrodynamic kick mechanism. 

The GW emission from aspherical flow in the Fryer et al. models sets in within 
~ 2-5 ms after bounce and, according to the authors, is mostly due to asymmetries 
caused by oscillations in the PNS core that are excited by the counteracting recoils 
due to aspherically ejected matter and neutrino momentum. However, and /ix of 
their models show no periodicity and appear uncorrelated. Considering the discussion 
on prompt convection in section 6, it appears not unlikely that at least part of the GW 
signal is due to prompt convection. The GW signals reach their maximum amplitudes 
of ~ 2-8 x 10~^^ (at 10 kpc) around 10-20 ms after bounce, decay thereafter and 
contain no large-amplitude GW memory. The absence of the latter is not discussed 
by Fryer et al., but may be due to the relatively small asphcricity of their explosions. 
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The authors provided neither GW spectra nor numbers for the energy emitted in 
GWs. 

9.3. Magnetic Fields 

As already pointed out in section 8, not only aspherical matter dynamics, but any kind 
of accelerated transport of energy, including magnetic-field energy, may lead to the 
emission of GWs. Kotake et al. [87] extended the Newtonian quadrupole formalism of 
wave generation to include the contribution from magnetic fields and introduced terms 
associated with Lorentz force and magnetic-field energy. They focussed on the GW 
signal of rotating core collapse and bounce and performed a limited set of Newtonian 
2D magneto-rotational core-collapse calculations with a finite-temperature nuclear 
EOS and a leakage scheme for neutrinos. They demonstrated that the magnetic 
contribution to the GW signal at core bounce remains smaller than ~ 1% of the 
matter contribution for cores with precollapse magnetization below ~ 10^^ G. The 
latter is already an extreme value, considering that garden-variety precollapse iron 
cores are likely to be weakly magnetized with B < 10^-10^°G [142]. Since matter and 
B-field are strongly coupled, the magnetic component of the GW signal exhibits time 
variations similar to those of the matter component. Interestingly, it leaves behind 
a secularly growing GW memory that, at 20 ms after bounce, dominates the overall 
signal with an amplitude of ^ 2 x 10~^^ (at 10 kpc) in Kotake et al.'s model with 
precollapse magnetic field of 5 x 10^^ G and central angular velocity fic,i = 4 rads^^. 
This GW memory may be due to the gradual build-up of the magnetic field energy 
in polar regions via an 17 dynamo (e.g., [17, 90, 113] and references therein), but was 
not studied in detail by Kotake et al. Also note that the models of Kotake et al. were 
not followed through explosion, hence do not exhibit a contribution to the GW signal 
and GW memory from aspherical outflows (see section 9.1). 

Obergaulinger et al. [90, 125] performed an extensive set of magneto-rotational 
core-collapse simulations with simplified microphysics, leading to prompt explosions 
with bipolar outflows (see section 9.1). They confirmed the previous result of 
Kotake et al. [87] that MHD effects on the dynamics and the magnetic contribution 
to the GW signal of core bounce remain small for cores with precollapse magnetic 
field strengths below ~ 10^^ G. For models with stronger precollapse magnetic fields, 
Obergaulinger et al. found significant magnetic contributions to the GW signal that 
grow after bounce due to the continuously increasing magnetic stresses in the polar 
regions and inside the PNS. The time variation in the magnetic component is much 
slower than in the matter component, peaking in frequency below ~ 100 Hz. In 
multiple models with strong magnetic fields (initial B > 10^^-10^^ G), the magnetic 
component eventually dominates the GW signal and leaves behind a GW memory that 
may have positive or negative sign, varying from model to model. The maximum GW 
signal amplitudes due to the magnetic component generally scale with the precollapse 
magnetic field strength. For their most extreme model, Obergaulinger et al. found a 
maximum magnetic GW signal amplitude of ^ 4.5 x 10"^^ (at 10 kpc) while most 
other models show one to two orders of magnitude smaller maximum amplitudes. Note 
that the GW memory due to the magnetic contribution linearly combines with the 
memory that is due to the bipolar outfiow of matter. Depending on the sign of the 
magnetic contribution, the total amplitude of the GW memory may be increased or 
decreased. 
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9.4- GWs from Collapse to a Black Hole 

Ordinary massive stars in the mass range from ^ 8 to ~ 100 Mq burn their nuclear 
fuel all the way to iron-group nuclei (or a mixture of oxygen, neon, and magnesium 
nuclei at the lowest masses) and their cores stay in hydrostatic quasi-equilibrium 
throughout their nuclear burning lives. The final iron or ONeMg cores are supported 
primarily by the pressure of degenerate electrons (Pe), secondarily by the thermal 
pressure of the heavy ions (Pion), and tertiarily by radiation pressure (Prad) such that 
Pe ^ Pion ^ Prad [4]. Hcncc, irou/ONcMg cores resemble white dwarfs and have 
effective Chandrasekhar masses in the range of 1.3 to 2.0 Mq [3, 137]. When 
such a core collapses, it separates into a subsonically contracting inner core and a 
supersonically infalling and consequently rarefying outer core. At bounce, the inner 
core has a mass of '--^ 0.5 to ^ 0.7 M© (depending on its average electron fraction Ye, 
average specific entropy s, and rotational configuration [108, 143]). The solid-core of 
the nuclear force and the resulting stiff EOS above nuclear density easily stabilizes 
the collapse of the inner core. Collapse always results in a PNS and direct collapse 
to a black hole (BH) without core bounce and a PNS phase occurs never in ordinary 
massive stars^. 

After bounce, outer core material accretes onto the PNS. The maximum baryonic 
NS mass+ that can be supported against gravity depends sensitively on the stiffness 
of the high-density EOS and ranges between ^ 1.5 Mq and ~ 2.5 — 3 Af© (e.g., 
[5, 228, 229]). Rapid differential rotation can enhance the maximum mass by up to 
^ 50% (e.g., [230]), while the thermal structure of the PNS has less influence on the 
maximum mass [229, 231]. 

In order to leave behind a neutron star, the SN explosion has to set in before the 
PNS has reached its maximum mass. Depending on the latter's actual value and on 
the accretion rate set by the structure of the progenitor, an explosion has to occur 
within ^ 1 — 2 s after bounce. If the explosion is weak and does not unbind the 
entire stellar envelope, fall-back accretion can still push the PNS over its mass limit 
at later times (e.g., [232]). When a PNS becomes gravitationally unstable before the 
explosion is launched, the SN engine is immediately shut off and the star ends its 
life as a coUapsar, possibly exploding in a gamma-ray burst if it possesses the needed 
amount and distribution of angular momentum [19, 20, 233]. 

Nonspherical PNS collapse to a BH results in a burst of GWs (1) due to the rapidly 
shrinking mass-quadrupole moment of the PNS during the collapse and (2) due to the 
quasi normal mode (QNM) ringing of the nascent BH that is initially distorted from 
its stationary Kerr shape and experiences subsequent distortion by accreting material. 

The GW emission from the collapse of a compact star to a BH was first considered 
in the late 1970's by Cunningham, Price, and Moncrief [234-236] whose work was 
later improved by Seidel and collaborators [237, 238] and more recently by Harada, 
Iguchi, and Shibata [239]. These authors employed a perturbative approach and 
studied the GW emission from nonspherical perturbations on a spherically-symmetric 
background. They found that the overall GW signal of BH formation is dominated by 

1 There is considerable confusion in the astrophysics community about this fact. This is in part 
due to the misleading use of the terms 'direct/prompt BH formation' and 'delayed BH formation' by 
Fryer, Woosley, and collaborators. In their context (e.g., [227]), these terms describe 'BH formation 
without SN explosion' and 'BH formation after SN explosion by fallback accretion', respectively. 
+ Astronomical observations generally measure the gravitational mass, the baryonic mass minus the 
mass-energy equivalent of the NS gravitational binding energy. 
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the QNM ringing of the formed BH and that most of the GWs are emitted from the 
fundamental quadrupole mode. If the perturbation is axisymmetric (i.e., primarily due 
to rotation) the emission is dominated by the £ = 2,m = 0,n ~ (200) multipolar 
component, while nonaxisymmetric perturbations show up as ^ = 2, m 7^ 0, n = 
modes. These modes emit at characteristic frequencies that depend only on the 
BH mass M and its dimensionless angular momentum parameter j = {cJ)/{GM'^). 
Berti et al. [240] provide fitting-formulae and coefficients that describe to better than 
~ 5% accuracy the BH QNM frequencies fimn- Focussing on the 200 (axisymmetric) 
and 220 (nonaxisymmetric, bar-like) QNMs, we find using their tables, 



Hence, a nonrotating BH with a mass of 2 Mq rings when perturbed in its 200 (or 220) 
mode with a frequency of ^ 6 kHz and this frequency decreases linearly as more matter 
is accreted. In the case of slow rotation (small j), one would expect a 'reverse chirp' 
from higher to lower frequencies in the GW signal as a BH is formed and accretes the 
massive stellar envelope in a failing core-collapse SN. With increasing spin, the QNM 
frequencies increase at different rates for different azimuthal mode number m and for 
rapid rotation, the decrease of the QNM frequencies by increasing BH mass could be 
partly compensated by increasing j. Initial and advanced LIGO-class detectors have 
their maximum sensitivity at frequencies below ~ 1 kHz, thus will require large signal 
amplitudes and emitted energies to detect even a galactic NS collapse event. 

By making estimates of the energy in the perturbations and its azimuthal m 
distribution, one can use equations 9 and 10 and analogous expressions for other m to 
estimate GW signal amplitudes. This can be done as outlined in [88] or by assuming 
that the GW signal is of sine-Gaussian shape as discussed, e.g., in [243]. However, this 
shall not be repeated here, since such estimates for spatial distribution and energy of 
perturbations cannot be made reliably. The latter highlights a fundamental problem of 
the perturbative approach: While the effects of a (small) perturbation on a background 
can be studied with great accuracy, the perturbation itself cannot be determined and 
must be put in by hand or provided as initial data from numerical simulations (e.g.. 



Direct numerical simulation of nonspherical BH formation using GR hydrody- 
namics coupled to spacetime curvature evolution allows for the self-consistent deter- 
mination of the GW signal. The first 2D simulations of BH formation were carried 
out in 1981 by Nakamura [245], but GWs could not be extracted due to numerical 
difficulties. The first waveforms and GW energy estimates from axisymmetric ro- 
tating stellar collapse to a BH were provided in 1985 by Stark and Piran [246] who 
used highly approximate initial data that were not in rotational equilibrium. More 
recently, Baiotti et al. [241, 242] presented the first GW signals from the 3D collapse 
of rigidly-rotating polytropic NS models in rotational equilibrium. They found that 
a collapsing rotating NS stays essentially axisymmetric and confirmed that the GW 
emission is primarily due to axisymmetric BH QNM oscillations. They also identi- 
fied a high-frequency, pre-BH formation component that is sensitive to the way the 
collapse of the polytropic NS is instigated. In figure 13, we present the GW signal 
(left panel) and GW power spectrum of the QNM ringing (right panel) extracted 




(10) 



(9) 



[244]). 
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Figure 13. Left: GW signal of tlic collapse of a rapidly rigidly rotating 
polytropic NS to a black hole in terms of the dominant even-parity metric 
perturbation as computed by Baiotti et al [241, 242]. Even-parity and odd- 
parity metric perturbations are related to the two physical GW polarizations 
according to h+ -ih^=^ T.i,m(QL " » /-oo <3L'i*')-2l'"", where ^ar"" is 
the —2 spin- weighted spherical harmonic [122, 241]. This plot corresponds to the 
bottom-right panel of figure 5 of [242] and is used by kind permission from the 
authors. Right: Power spectrum of the BH ring-down part of the GW signal for 
the same model the units arc chosen in such a way that the peak of the power 
spectrum is approximately 1. The dotted line marks the BH QNM frequency 
obtained by perturbation theory. This plot corresponds to the bottom-right panel 
of figure 6 of [242] and is used by kind permission from the authors. 

from the rapidly rigidly rotating collapse model D4 with M = 1.86 and j = 0.54 
of Baiotti et al. [242]. While the early part of the waveform is dominated by de- 
tails of the matter dynamics and may be sensitive to the EOS and the NS angular 
momentum distribution, the BH QNIVI ringing depends only on the mass M and an- 
gular momentum parameter j, hence, will be qualitatively and quantitatively similar 
in simulations that take into account a more realistic matter treatment and differen- 
tial rotation. For their model suite, Baiotti et al. [241, 242] give a maximum energy 
emission of i?GW = 1-45 x 10^^{M/Mq)Mqc'^ which obtains in their most rapidly 
rotating model D4 and scales roughly as at small j. Convolving their GW signals 
with detector noise curves as discussed in [34], they find an upper-limit characteristic 
strain he of ~ 5.5 x 1O'^^(M/M0) at 10 kpc and at a characteristic frequency /c of 
^ 530 Hz for initial LIGO detectors. These numbers suggest that BH formation in a 
galactic core-collapse event with rapid rotation may be marginally detectable already 
with initial LIGO-class detectors. 

Simulations that track BH formation in a microphysically detailed multi-D model 
of a failing core-collapse SN remain to be carried out. Such simulations are necessary 
for the self-consistent study of the early, pre-QNlVI ringing GW signal of PNS collapse 
as well as for the GW signal of BH QNIVI oscillations excited by long-term accretion 
(but see [247] for a perturbative treatment of the latter). 
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Table 4. List of optically discovered core-collapse SNe that occurred 
within 5 Mpc from Earth between January 1, 2002 and August 31, 
2008. The table is based on the comprehensive online SN listing at 
http://www.cfa.harvard.edu/iau/lists/Sup6rnovae.htinl. Given are the SN 
name, the name of the host galaxy, the date of discovery, the core-collapse SN 
spectral subtype (see, e.g., [250]), and the approximate distance. Note that the 
date of astronomical discovery always postdates the actual explosion. 



SN 


Host Galaxy 


Date 


Type 


Distance 


2008bk 


NGC 7793 


20080325 [251] 


II-P 


~ 3.9 [252] 


2005af 


NGC 4945 


20050208 [253] 


II-P 


~ 3.6 [252] 


2004dj 


NGC 2403 


20040731 [254] 


II-P 


~ 3.3 [252] 


2004am 


M 82 


20040305 [255] 


II-P 


~ 3.5 [256] 


2002kg 


NGC 2403 


20021026 [257] 


Iln 


~ 3.3 [252] 



10. Nearby Core-Collapse SNe, the recent SN 2008bk, and their 
Detectability in Gravitational Waves 

Optimistic estimates of the core-collapse SN rate in the IVIilky Way and the close-by 
Small and Large ]VIagellanic Clouds predict one core-collapse SN in ^ 30 — 50 years and 
even for the entire local group of galaxies, including M31 at 0.8 IVIpc, one core-collapse 
SN in ~ 20 years is an optimistic estimate (see, e.g., [37] and the compilation of rate 
estimates and references in [42]). This rate stays roughly constant until a distance 
of ^ 3Mpc. The galaxies of the M81 group and neighboring groups with high star- 
formation rates that are located at 3-5 IVIpc from Earth increase the core-collapse SN 
rate to an optimistic 1 core-collapse SN in ^ 2 years within ~ 5 IVIpc [38, 248]. The 
next significant increase in the SN rate occurs when the outskirts of the Virgo cluster 
are beginning to contribute at 7-10 Mpc [38]. 

Since LIGO/GEO600 [26, 27] science operations began in 2002, five core-collapse 
SNe have been discovered optically within 5 JVIpc from Earth, a number roughly 
consistent with the above quoted rate estimate. Unfortunately, all events occurred 
outside of LIGO/GEO600 and VIRGO [28] science runs. However, the very recent 
(as of August 2008) nearby core-collapse SN, SN 2008bk, occurred while GEO600 and 
the 2-km LIGO Hanford interferometer (H2) were taking data in "astrowatch" mode 
[249]. All other LIGO-class detectors, including VIRGO, were offiine for upgrades. 

SN 2008bk was found on IVIarch 25, 2008 in NGC 7793 [251], a sphal galaxy 
approximately 3.9 IVIpc away from Earth. Based on the lightcurve evolution of similar 
type II-P ("plateau," see, e.g., [250]), 2008bk exploded ~ 20-36 days before its 
discovery [251]. The progenitor star appears to have been on the low-mass end of 
the massive-star population and may have had a mass between 8.5 ± 1 Mq [258, 259]. 

Since LIGO H2 and GEO600 were most hkely taking data when SN 2008bk 
exploded, it is interesting to estimate the detectability by H2 and GEO600 for the 
GW emission mechanisms discussed in in this review. Given the high rate of core- 
collapse SNe between ~ 3 and ~ 5 IVIpc it is also useful to consider their detectability 
by other current and future GW observatories. 

Since waveforms for some of the GW emission mechanisms are available, we 
compute single-detector optimal (assuming perfect orientation) matched-filtering 
signal-to-noise ratios (SNR), using 
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Table 5. Upper-limit SNR estimates for SN 2008bk at 3.9 Mpc using theoretical 
noise power spectral densities for advanced LIGO 4-km interferometers (LIGO 2 
4km), LIGO 1 4-km interferometers (LIGO Ll/Hl), VIRGO, and LIGO S5 data 
for the H2 and GEO600 interferometers. Considered are various representative 
models for GW emission by rotating core collapse and bounce, nonaxisymmetric 
rotational instabilities, and PNS core g-modes. Behind each model name a paper 
reference is provided. For some models, ^^-scaled SNRs are given as estimates 
for longer GW emission than tracked by the model calculations. This is denoted 
by (xn) behind the model name. 



Process 


Model 


LIGO 2 


LIGO Ll/Hl 


LIGO H2 


GEO600 


VIRGO 






4 km 


4 km 


2 km 


600 m 


3 km 


Rotating Collapse 


S11A2013 


0.124 


0.008 


0.005 


0.001 


0.009 


&c Bounce 


S20A2O09 


0.130 


0.008 


0.006 


< 0.001 


0.010 


[108] 


S40A3O12 


0.214 


0.024 


0.013 


< 0.001 


0.018 


Rotational Instability 


S20A2B4 


0.319 


0.021 


0.014 


0.003 


0.022 


[42, 115, 119] 


S20A2B4 {x5) 


0.713 


0.047 


0.031 


0.007 


0.049 


PNS g-modes 


sll.2 


0.147 


0.006 


0.005 


0.002 


0.009 


[22, 23] 


sl5.0 


0.454 


0.021 


0.015 


0.006 


0.027 


and section 7.1 


s25.0 


0.612 


0.029 


0.020 


0.007 


0.037 




s25.0 (x2) 


0.866 


0.041 


0.029 


0.009 


0.052 




s25WW 


5.331 


0.217 


0.151 


0.057 


0.328 



where S{f) is the one-sided detector noise power spectral density in units of Hz ' 
and h is the Fourier transform of the wave signal, computed via [34, 191] 

/oo 
e^^'f*h+^yit)dt . (12) 
-oo 

Using the characteristic strain /ichar(/) defined by equation 2, equation 11 can be 
rewritten [191] to 

f°° h"^ ( f] 

(SNR)2p,,„,i = / dln/-fM^, (13) 

where /irms(/) is the dimensionless detector rms noise given by /irms(/) = \/ fS{f). 
Note that the expression for the optimal SNR given here is a factor of 3/2 smaller 
than the expression used in [108] and derived in [34]. See [191] for details. 

For computing SNRoptimah we employ noise power spectral densities S{f) for H2 
and GEO600 from the LIGO/GEO600 S5 runs. The H2 S{f) data have so far not been 
published and were kindly made available to us by M. Landry [260]. The GEO600 
S{f) data were obtained from the GEO600 website [27]. In addition, we compute 
SNRoptimai based on the design S{f) for VIRGO [28], LIGO Livingston/Hanford 1 
[190], and advanced LIGO 4-km interferometers (LIGO 2) in burst mode [190]. 

In table 5, we list SNRoptimai for all considered detectors and for models focussing 
on GW signals from rotating collapse and bounce, nonaxisymmetric rotational 
instabilities, and PNS 5-modes at the distance of SN 2008bk. GW emission from 
convectivc overturn and the other emission mechanisms considered in this review are 
not included in table 5, since their SNRoptimai would be of order 10~^ and below at a 
distance of ~ 4 Mpc even for advanced LIGOs (see, e.g., [97] and sections 6 and 9). 

The first model given for each emission process in table 5 is an "average emitter" 
and the subsequent models range from "strong" to "extreme" . Since in many 
computed models the simulation was stopped before the GW emission had subsided, 
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we also provide y/n-scaled (where n is the factor by which the emission is prolonged) 
SNRoptimai for a subset of models. 

Based on the SNRoptimai listed in table 5, we surmise that LIGO H2 and GEO600 
had no chance of seeing GWs emitted in SN 2008bk. Even the most extreme models 
yield optimal SNRs below 0.2 at a source distance of 3.9 Mpc. For a detection, a 
SNR significantly greater than 1 and probably in the range of 5 — 7 would be needed 
(e.g., [97, 191]). Once in operation and at design sensitivity, advanced LIGO could 
constrain the GW emission from PNS g-modes out to the distance of SN 2008bk and 
beyond. All other emission processes require a closer core-collapse event (or more 
sensitive future detectors) to be observable. 

11. Putting things together: Summary and a Conjecture 

The GW signature of core-collapse SNe is rich and multi-faceted. The aim of this 
review was to provide an overview and summary of the current knowledge about the 
various GW emission processes that may occur in a core-collapse event. 

We have outlined and discussed rotating core collapse and bounce, nonaxisym- 
metric rotational instabilities, postbounce convective overturn, and non-radial PNS 
pulsations as the prominent candidate processes whose multi-D dynamics are likely to 
emit GWs. In addition, we have summarized the emission characteristics of the GW 
signals associated with anisotropic neutrino emission, aspherical outflows, magnetic 
stresses, global precollapse asymmetries, and the collapse of hypermassive PNSs to 
BHs. 

All of the listed processes are burst emitters of GWs, but range in temporal and 
spectral characteristics from short, fcw-ms bursts with rather narrow-band emission 
(e.g., rotating collapse and bounce) to broadband emission lasting, perhaps, for 
seconds (e.g., convection, in particular PNS convection). In addition, a number of 
the emission processes produce bursts with GW memory, leaving behind essentially 
zero-frequency, permanent distortions of spacetime. 

Rotating core collapse and bounce is the most extensively investigated and 
arguably the best understood source of GWs in the core collapse context. 2D and 
3D simulations in conformally-flat and full GR were carried out and included the 
dominant microphysical aspects, results of different approaches agree, and waveform 
catalogs arc available [126, 155] to the GW data analysis community which is beginning 
to employ them (e.g., [261, 262]). Nevertheless, uncertainties remain. Without 
doubt, the biggest is the lack of presupernova models from multi-D stellar evolution 
simulations. Present simulations have to rely upon profiles from ID stellar evolution 
codes that take into account crucial multi-D effects such as convection and rotation 
in only highly approximate ways. Core collapse investigations try to deal with this 
by large-scale parameter studies, but until multi-D stellar evolution models become 
available and reliable, even extensive parameter studies could, in principle, be missing 
the real precollapse configurations that obtain in nature. 

In order to produce a burst of GWs at core bounce comfortably detectable by 
advanced LIGOs throughout the Milky Way, a precollapse iron core needs ([108] and 
section 4) to be spinning very rapidly with central periods -Po 6 s (corresponding 
to central angular velocities > 1 rads~^). Initial LIGOs might see cores with 
-fo ^ 3 s. Perhaps only '--^ 1% or less of the massive star population rotates with 
such short periods at the onset of core collapse [144], while the vast majority is likely 
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spinning much more slowly with periods of 30— 100 s as currently predicted [141, 142]. 

Rapidly spinning iron cores with the above stated short central periods produce 
millisecond-period PNSs with rotational energy of order 10 B of which a fraction 
would be sufficient to power an energetic SN explosion. Moreover, such PNSs have 
strongly differentially-rotating outer cores (e.g., [108, 141, 263]) and are likely to 
become subject to nonaxisymmetric rotational instabilities which would lead to the 
GW emission characteristics outlined in section 5. As discussed by a number of groups 
(e.g., [16, 17, 20, 90, 117, 165]), magnetic field amplification can draw from rotational 
energy. In particular, Burrows et al. [17] found that cores with Pq ^ 2 — 4 s develop 
MHD-driven jet-like SN explosions*. Hence, there is a close link between a detectable 
GW burst from core bounce and MHD-driven explosions. The same may be true for 
3D rotational instabilities and their GW signal, but it is presently not known in which 
way 3D hydrodynamic and MHD instabilities interact in a postbounce SN core. 

Convection and the SASI are very likely to be present in the postbounce pre- 
explosion phase of the vast majority of core-collapse events. Rapid rotation damps 
convection and a very early onset of explosion, probably relevant only in the case of 
ONeMg cores in the lowest-mass massive stars [2, 12], can shut off convection before 
it is fully developed. 

In contrast to the GW signal of rotating collapse and bounce, the GW emission 
characteristics of postbounce convective overturn and the SASI have been investigated 
by only a few studies and the systematics of the GW signal with variations in 
progenitor structure, rotational configuration etc. have not been investigated in detail. 
Furthermore, owing to the stochastic nature of both large-eddy and small-scale 
turbulent convection, GW templates that could be expected to match real signals 
cannot be predicted. It will most likely take an advanced LIGO-class detector to detect 
for a galactic event the broadband ^ 0.5—1.5 s GW signal emitted by convection/SASI 
before explosions and the longer-duration, but lower-amplitude GW signal from PNS 
convection. 

The explosion scenario currently favored for slowly rotating and nonrotating 
massive stars - the neutrino mechanism (e.g., [2, 4, 7]) - relies on the deposition of 
sufficient energy by neutrinos in the postshock region to revive the stalled shock and 
unbind the stellar envelope. A multitude of multi-D simulations have demonstrated 
that convection and SASI are key ingredients for neutrino-driven core-collapse SN 
explosions (e.g., [9, 10] and references therein). The acoustic mechanism, recently 
proposed by Burrows et al. [21, 22] and intensely debated in the community [9, 24, 25], 
is an alternative to the neutrino mechanism and relies on the excitation by accretion 
and turbulence of large-amplitude PNS core g-modes and their damping by strong 
sound waves that deposit energy in the postshock region, eventually leading to 
explosion. The GW signature of an acoustically-driven core-collapse SN explosion 
would be dominated by the GWs emitted from the quadrupole components of the 
PNS core pulsations. Already first-generation LIGOs should be able to see such a 
strong signal throughout the Milky Way. Advanced detectors may be able to put 
significant constraints on the signal strength out to 3 — 5 Mpc where the SN rate is 
favorably high. 

In table 6 we summarize the prominent GW emission processes that are expected 
to be active in a core-collapse event and contextualize them with the three core-collapse 

* Provided either that the magneto-rotational instability (MRI) [164] occurs as envisioned or that 
the precollapse core magnetic field is at least of the order of 10^^ G. 
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Table 6. Overview on prominent GW emission processes in core-collapse SNe and 
their possible emission 'levels' in the context of the magneto-rotational (MHD; 
e.g. [17]), the neutrino (e.g., [9]), and the acoustic SN mechanism (e.g., [21, 22]). 
For a galactic SN, 'strong' corresponds to 'probably detectable by initial and 
advanced LIGO', 'weak' means 'probably marginally detectable by advanced 
LIGO', and 'none' means 'absent or probably not detectable by advanced LIGO'. 
The three considered explosion mechanisms are likely to have mutually exclusive 
GW signatures and could be distinguished by the mere detection or non-detection 
of a GW signal without the solution of the full inverse problem. Note that the GW 
signal due to anisotropic neutrino emission, though present in all three mechanism, 
is not considered, since its low-frequency character severely limits its detectability. 



GW Emission 
Process 


Pc 

MHD Mechanism 
(rapid rotation) 


tential Explosion Mecli 
Neutrino Mechanism 
(slow/no rotation) 


mism 
Acoustic Mechanism 
(slow/no rotation) 


Rotating Collapse 
and Bounce 


strong 


none / weak 


none / weak 


3D Rotational 
Instabilities 


strong 


none 


none 


Convection 
& SASI 


none / weak 


weak 


weak 


PNS g-modes 


none / weak 


none/ weak 


strong ^ 



SN explosion mechanisms by semi-quantitatively assigning 'emission strengths' for a 
galactic event {D « 10 kpc) based on the GW emission estimates collected in this 
review. In the scale we assume, 'none' refers to 'absent or probably not detectable 
by advanced LIGOs (SNR <C 1)', 'weak' means 'probably only marginally detectable 
by advanced LIGOs (SNR > 1),' 'strong' reflects 'probably detectable by initial and 
advanced LIGOs (SNR > 5 — 7).' The roughness of the scale is intended to emphasize 
the independence of the overall argument from variations in quantitative details. 

The MHD mechanism is limited to rapidly rotating cores, hence will involve strong 
GW emission from core bounce and, perhaps, nonaxisymmetric rotational instabilities. 
The acoustic mechanism, on the other hand, may work best in nonrotating or slowly 
rotating cores and will emit a strong GW signal from PNS pulsations. The neutrino 
mechanism, also probably most relevant in slowly or nonrotating cores, will very likely 
have convection/SASI as its strongest source of GWs. 

Based on the above discussion and on table 6, we conjecture that the GW 
signatures of the neutrino, MHD, and acoustic core-collapse SN mechanisms are 
mutually exclusive. Hence, for a galactic SN and even initial LIGOs, the mere detection 
of a GW signal and its association with an emission process and, in fact, also the non- 
detection of GWs, have the potential of strongly constraining the way massive stars 
explode. 

Unfortunately, galactic core-coUapsc SNc arc quite rare events, occurring at a 
rate of one in a few decades. Initial LIGO-class detectors cannot see core-collapse 
SNe outside the Milky Way. Advanced detectors, at least as currently planned [190], 
may be able to see most core-collapse SNe throughout the local group of galaxies 
{D « 1 Mpc). This, however, would increase the observable event rate by not more 
than a factor of 2 [37, 38]. Third-generation GW observatories (e.g. the envisioned 
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EURO detector [264] ) reaehing out to at least 3 — 5 Mpc will be neeessary for detailed 
GW astronomy of eore-coUapse SNe to beeome possible. 
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